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Quantum state engineering, i.e., active control over the coherent 
dynamics of suitable quantum mechanical systems, has become a 
fascinating perspective of modern physics. With concepts devel- 
oped in atomic and molecular physics and in the context of NMR, 
the field has been stimulated further by the perspectives of quan- 
tum computation and communication. For this purpose a number 
of individual two-state quantum systems (qubits) should be ad- 
dressed and coupled in a controlled way. Several physical realiza- 
tions of qubits have been considered, inch trapped ions, NMR, and 
quantum optical systems. For potential applications such as logic 
operations, nano-electronic devices appear particularly promising 
because they can be embedded in electronic circuits and scaled up 
to large numbers of qubits. 

Here we review the quantum properties of low-capacitance 
Josephson junction devices. The relevant quantum degrees of free- 
dom are either Cooper pair charges on small islands or fluxes in 
ring geometries, both in the vicinity of degeneracy points. The co- 
herence of the superconducting state is exploited to achieve long 
phase coherence times. Single- and two-qubit quantum manipu- 
lations can be controlled by gate voltages or magnetic fields, by 
methods established for single-charge systems or the SQUID tech- 
nology. Several of the interesting single-qubit properties, incl. co- 
herent oscillations have been demonstrated in recent experiments, 
thus displaying in a spectacular way the laws of quantum mechanics 
in solid state devices. Further experiments, such as entanglement 
of qubit states, which are crucial steps towards a realization of logic 
elements, should be within reach. 

In addition to the manipulation of qubits the resulting quan- 
tum state has to be read out. For a Josephson charge qubit this 
can be accomplished by coupling it capacitively to a single-electron 
transistor (SET). To describe the quantum measurement process 
we analyze the time evolution of the density matrix of the coupled 
system. As long as the transport voltage is turned off, the transis- 
tor has only a weak infiuence on the qubit. When the voltage is 
switched on, the dissipative current through the SET destroys the 
phase coherence of the qubit within a short time. The measurement 
is accomplished only after a longer time, when the macroscopic sig- 
nal, i.e., the current through the SET, resolves different quantum 
states. At still longer times the measurement process itself destroys 
the information about the initial state. Similar scenarios are found 
when the quantum state of a fiux qubit is measured by a dc-SQUID, 
coupled to it inductively. 
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I. INTRODUCTION 

The interest in 'macroscopic' quantum effects in low- 
capacitance Josephson junction circuits has persisted for 
many years. One of the motivations was to test whether 
the laws of quantum mechanics apply in macroscopic 
systems, in a Hilbert space spanned by macroscopically 
distinct states (Leggett 1987). The degrees of freedom 
studied were the phase difference of the superconduct- 
ing order parameter across a junction, or the flux in a 
superconducting ring (SQUID) geometry. Various quan- 
tum phenomena, such as macroscopic quantum tunneling 
(MQT) and resonance tunneling were demonstrated (see 
e.g. Voss & Webb 1981, Martinis et al. 1987, Rouse et al. 
1995). On the other hand, despite experimental efforts 
(e.g. Tesche 1990) coherent oscillations of the flux be- 
tween two macroscopically distinct states (macroscopic 
quantum coherence, MQC) had not been observed. 

The field received new attention recently, after it was 
recognized that suitable Josephson devices may serve as 
quantum bits (qubits) in quantum information devices, 
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and that quantum logic operations can be performed by 
controlling gate voltages or magnetic fields (e.g. Bouch- 
iat 1997, Shnirman et al. 1997, Averin 1998, Makhlin 
et al. 1999, Nakamura et al. 1999, Mooij et al. 1999, loffe 
et al. 1999). In this context, as well as for other con- 
ceivable applications of quantum state engineering, the 
experimental milestones are the observation of quantum 
superpositions of macroscopically distinct states, of co- 
herent oscillations, and of entangled quantum states of 
several qubits. For Josephson devices first successful ex- 
periments have been performed. These systems can be 
fabricated by established lithographic methods, and the 
control and measurement techniques are quite advanced. 
They further exploit the coherence of the superconduct- 
ing state, which helps achieving sufficiently long phase 
coherence times. 

Two alternative realizations of quantum bits have been 
proposed, based on either charge or phase (flux) degrees 
of freedom. In the former the charge in low-capacitance 
Josephson junctions is used as quantum degree of free- 
dom, with basis states differing by the number of Cooper 
pair charges on an island. These devices combine the co- 
herence of Cooper pair tunneling with the control mecha- 
nisms developed for single-charge systems and Coulomb- 
blockade phenomena. The manipulations can be ac- 
complished by switching gate voltages (Shnirman et al. 
1997); designs with controlled inter-qubit couplings were 
proposed (Averin 1998, Makhhn et al. 1999). Exper- 
imentally, the coherent tunneling of Cooper pairs and 
the related properties of quantum mechanical superposi- 
tions of charge states has been demonstrated (Bouchiat 
1997, Nakamura et al. 1997). Most spectacular are recent 
experiments of Nakamura et al. (1999), where the quan- 
tum coherent oscillations of a Josephson charge qubit 
prepared in a superposition of eigenstates were observed 
in the time domain. We describe these systems, concepts 
and results in Section ||. 

The alternative realization is based on the phase of a 
Josephson junction or the flux in a ring geometry near 
a degeneracy point as quantum degree of freedom (see 
e.g. Mooij et al. 1999, loffe et al. 1999). In addition to 
the earlier experiments, where macroscopic quantum tun- 
neling had been observed (Voss & Webb 1981, Martinis 
et al. 1987, Rouse et al. 1995), the groups in Delft and 
Stony Brook (van der Wal et al. 2000, Friedman et al. 
2000) demonstrated recently by spectroscopic measure- 
ments the flux qubit's eigenenergies, they observed eigen- 
states which are superpositions of different flux states. 



and new efforts are made to observe the coherent oscil- 
lation of the flux between degenerate states (Mooij et al. 
1999, Friedman et al. 2000, CosmeUi et al. 1998). We 
will discuss the quantum properties of flux qubits in Sec- 
tion 



III 



To make use of the quantum coherent time evolution 
it is crucial to find systems with intrinsically long phase 
coherence times and to minimize external sources of de- 
phasing. The latter can never be avoided completely 
since, in order to perform the necessary manipulations, 
one has to couple to the qubits, e.g., by attaching exter- 
nal leads. Along the same channels as the signal (e.g., 
gate voltages) also the noise enters the system. However, 
by operating at low temperatures and choosing suitable 
coupling parameters, these dephasing effects can be kept 
at an acceptable level. We provide estimates of the phase 
coherence time in Section [V. 



^ Since computational applications are widely discussed, we 
employ here and below frequently the terminology of quantum 
information theory, referring to a two-state quantum system 
as qubit and denoting unitary manipulations of its quantum 
state as quantum logic operations or gates. 



In addition to controlled manipulations of qubits, 
quantum measurement processes are needed, e.g., to read 
out the final state of the system. In our quantum mechan- 
ics courses we learned to express the measurement pro- 
cess as a "wave function collapse", i.e., as a non-unitary 
projection, which reduces the quantum state of the qubit 
to one of the possible eigenstates of the observed quantity 
with state-dependent probabilities. However, in reality 
any measurement is performed by a device which itself 
is realized by a physical system, suitably coupled to the 
measured quantum system and with a macroscopic read- 
out variable. Its presence, in general, disturbs the quan- 
tum manipulations. Therefore the dissipative processes 
which accompany the measurement should be switched 
on only when needed. 

An example is provided by a normal-state single- 
electron transistor (SET) coupled capacitively to a single- 
Cooper pair box. This system is widely used as an 
electro-meter in classical single-charge systems. We de- 
scribe in Section how a SET can also be used to read 
out the quantum state of a charge qubit. For this purpose 
we study the time evolution of the density matrix of the 
coupled system (Shnirman & Schon 1998). During quan- 
tum manipulations of the qubit the transport voltage of 
the SET is turned off, in which case it acts only as an ex- 
tra capacitor. To perform the measurement the transport 
voltage is turned on. In this stage the dissipative current 
through the transistor dephases the the state of the qubit 
rapidly. This current also provides the macroscopic read- 
out signal for the quantum state of the qubit. However, 
it requires a longer 'measurement time' until the noisy 
signal resolves different qubit states. Finally, on the still 
longer 'mixing time' scale, the measurement process it- 
self destroys the information about the initial quantum 
state. 

Many results and observations made in the context 
of the normal state single-electron transistor apply also 
to other physical systems, e.g., a superconducting SET 
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(SSET) coupled to a charge qubit (Averin 20006, Cot- 
tet et al. 2000) or a dc-SQUID monitoring as a quan- 
tum magneto-meter the state of a flux qubit (see e.g. 
Mooij et al 1999, Friedman et al. 2000, Averin 20006). 
The results can also be compared to the nonequilib- 
rium dephasing processes discussed theoretically (Levin- 
son 1997, Aleiner et al. 1997, Gurvitz 1997) and demon- 
strated experimentally by Buks et al. (1998). 

One of the motivations for quantum state engineer- 
ing with Josephson devices is their potential application 
as logic devices and quantum computing. By exploit- 
ing the massive parallelism of the coherent evolution of 
superpositions of states quantum computers could per- 
form certain tasks which no classical computer can do in 
acceptable times (Bennett 1995, Barenco 1996, DiVin- 
cenzo 1995, Aharonov 1998). In contrast to the develop- 
ment of physical realizations of qubits and gates, i.e., the 
"hardware", the theoretical concepts of quantum com- 
puting, the "software" , are already rather advanced. As 
an introduction, and in order to clearly define the goals, 
we present in Appendix ^ an ideal model Hamiltonian 
with sufficient control to perform all the needed manip- 
ulations. (We can mention that the Josephson junction 
devices come rather close to this ideal model.) Then in 
Appendix |b| we show by a few representative examples 
how these manipulations can be combined for useful com- 
putations. 

Various other physical systems have been suggested 
as possible realizations of qubits and gates. We men- 
tion ions in electro-magnetic traps manipulated by laser 
irradiation (Cirac & ZoUer 1995, Monroe et al. 1995), 
NMR on ensembles of molecules in liquids (Gershenfeld 
& Chuang 1997, Cory et al. 1997), and cavity QED sys- 
tems (Turchette et al. 1995). In comparison, the men- 
tioned Josephson systems are more easily embedded in 
electronic circuits and scaled up to large registers. Ultra- 
small quantum dots with discrete levels, and in particu- 
lar, spin degrees of freedom embedded in nano-structured 
materials are candidates as well. They can be manipu- 
lated by tuning potentials and barriers (Loss & DiVin- 
cenzo 1998, Kane 1998). InAppendix|| we describe these 
alternative solid state realizations of qubits and look at 
their advantages and drawbacks. Because of the difhcul- 
ties of controlled fabrication their experimental realiza- 
tion is still at a very early stage. 



II. JOSEPHSON CHARGE QUBIT 

A. Superconducting charge box as a quantum bit 

In this section we describe the properties of low- 
capacitance Josephson junctions, where the charging en- 
ergy dominates over the Josephson coupling energy, and 
discuss how they can be manipulated in a quantum co- 
herent fashion. Under suitable conditions they provide 



physical realizations of qubits with two states differing by 
one Cooper pair charge on a small island. The necessary 
one-bit and two-bit gates can be performed by control- 
ling applied gate voltages and magnetic fields. Different 
designs will be presented which differ in complexity, but 
also in the accuracy and flexibility of the manipulations. 



n 
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FIG. 1. A Josephson charge qubit in its simplest design 
formed by a superconducting single-charge box. 

The simplest Josephson junction qubit is shown in 
Fig. 1^. It consists of a small superconducting island 
( "box" ) with n excess Cooper pair charges (relative to 
some neutral reference state), connected by a tunnel junc- 
tion with capacitance Cj and Josephson coupling energy 
Ej to a superconducting electrode. A control gate volt- 
age Vg is coupled to the system via a gate capacitor Cg. 
Suitable values of the junction capacitance, which can be 
fabricated routinely by present-day technologies, are in 
the range of Femtofarad and below, Cj < 10~^^ F, while 
the gate capacitances can be chosen still smaller. The 
relevant energy scale, the single-electron charging energy 
Ec = e^/2(Cg -f Cj), which depends on the total canac- 
itance of the island, is then in the range of Kelvin or 
above, Eq > 1 K. The Josephson coupling energy i?j is 
proportional to the critical current of the Josephson junc- 
tion (see e.g. Tinkham 1996). Typical values considered 
here are in the range of 100 mK. 

We choose a material such that the superconducting 
energy gap A is the largest energy in the problem, larger 
even than the single-electron charging energy. In this 
case quasiparticle tunneling is suppressed at low temper- 
atures, and a situation can be reached where no quasi- 
particle excitation is found on the island El. Under these 



^Throughout this review we frequently use temperature 
units for energies. 

■^In the ground state the superconducting state is totally 
paired, which requires an even number of electrons on the 
island. A state with an odd number of electrons necessar- 
ily costs an extra quasi-particle energy A and is exponen- 



3 



conditions only Cooper pairs tunnel - coherently - in the 
superconducting junction, and the system is described by 
the Hamiltonian: 



— E-i cos 9 



(2.1) 



Here, n is the number operator of (excess) Cooper pair 
charges on the island, and G, the phase of the supercon- 
ducting order parameter of the island, is its quantum me- 
chanically conjugate, n = —ihd/d{hQ). The dimension- 
less gate charge, Ug = CgV^/2e, accounts for the effect of 
the gate voltage and acts as a control parameter. Here 
we consider systems where the charging energy is much 
larger than the Josephson coupling energy, Eq ^ Ej. In 
this regime a convenient basis is formed by the charge 
states, parameterized by the number of Cooper pairs n 
on the island. In this basis the Hamiltonian (2T) reads 



?i = ^{4i?c(n-ng)2|n)(n| 

n 

-ii?j(|n)(n+l| + |n + l)(n|)} . (2.2) 




FIG. 2. The charging energy of the superconducting electron 
box is shown as a function of the gate charge rig for differ- 
ent numbers of extra Cooper pairs n on the island (dashed 
parabolas). Near degeneracy points the weaker Josephson 
coupling mixes the charge states and modifies the energy of 
the eigenstates (solid lines). In the vicinity of these points the 
system effectively reduces to a two-state quantum system. 



For most values of Ug the energy levels are dominated 
by the charging part of the Hamiltonian. However, when 
Ug is approximately half-integer and the charging ener- 
gies of two adjacent states are close to each other (e.g., 
at Vg = Vdcg = e/Cg), the Josephson tunneling mixes 
them strongly (see Fig. 0). We concentrate on such a 



tially suppressed at low T. This 'parity effect' has been 
established in experiments below a crossover temperature 
T* »i A/(A;b In A'^off), where A^otr is the number of electrons in 
the system near the Fermi energy (Tuominen et al. 1992, La- 
farge et al. 1993, Schon & Zaikin 1994, Tinkham 1996). For 
a small island, T* is typically one order of magnitude lower 
than the superconducting transition temperature. 



voltage range near a degeneracy point where only two 
charge states, say n = and n = 1, play a role, while 
all other charge states, having a much higher energy, can 
be ignored. In this case the superconducting charge box 
( ^.iD reduces to a two-state quantum system (qubit) with 
Hamiltonian which can be written in spin-^ notation as 



1 



1 



Ctrl 



(2.3) 



The charge states n — and n = 1 correspond to the spin 
basis states |t) = (J) and ||) = (5), respectively. The 
charging energy splitting, which is controlled by the gate 
voltage, corresponds in spin notation to the z-component 
of the magnetic field 



B, = SE,i, = AEc{l-2ng), 



(2.4) 



while the Josephson energy provides the x-component of 
the effective magnetic field 



Bx = E. 



J • 



(2.5) 



For later convenience we rewrite the Hamiltonian as 

Tictri = -^AE{ri) (cosjy + sin?? ax) , (2.6) 

where the mixing angle i] = ta.ii~^{Bx/ B^) determines 
the direction of the effective magnetic field in the x-z- 
plane, and the energy splitting between the eigenstates 
is AE{ri) = ■\/B|"+~Bf = Ej/ sin 77. At the degeneracy 
point, 77 = 7r/2, it reduces to Ej. The eigenstates are 
denoted in the following as |0) and |1). They depend on 
the gate charge Ug as 



|0) = cos I IT) + sin I U) 
|l>=-sin||T>+cos||i) 



(2.7) 



We can further express the Hamiltonian in the basis 
of eigenstates. To avoid confusion we introduce a second 
set of Pauli matrices, p, which operate in the basis |0) 
and |1), while reserving the operators cr for the basis of 
charge states |t) and ||). By definition the Hamiltonian 
then becomes 



H = -^AEi7^)p, 



(2.8) 



The Hamil toni an (2.3) is similar to the ideal single- 
qubit model ( |Al|) presented in Appendix Ideally the 
bias energy (effective magnetic field in z-direction) and 
the tunneling amplitude (field in x-direction) are con- 
trollable. However, at this stage we can control - by the 
gate voltage - only the bias energy, while the tunneling 
amplitude has a constant value set by the Josephson en- 
ergy. Nevertheless, by switching the gate voltage we can 
perform the required one-bit operations (Shnirman et al. 
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1997). If, for example, one chooses the idle state far to 
the left from the degeneracy point, the eigenstates |0) and 
|1) are close to |t) and ||), respectively. Then, switching 
the system suddenly to the degeneracy point for a time 
At and back produces a rotation in spin space, 



C^i-bit(a) = exp 



cos § i sin ^ 
i sin ^ cos ^ 



(2.9) 



where a = EjAt/h. Depending on the value of At, a spin 
flip can be produced, or, starting from |0), a superposi- 
tion of states with any chosen weights can be reached. 
(This is exactly the operation performed in the experi- 
ments of Nakamura et al. (1999); see Subsection [I.D| ). 
Similarly, a phase shift between the two logical states can 
be achieved by changing the gate voltage rig for some time 
by a small amount, which modifies the energy difference 
between the ground and excited states. 
Several remarks are in order: 

(1) Unitary rotations by and are sufficient for all 
manipulations of a single qubit. By using a sequence 
of no more than three such elementary rotations we can 
achieve any unitary transformation of a qubit's state. 

(2) The example presented above, with control of B^ only, 
provides an approximate spin flip for the situation where 
the idle point is far from degeneracy and Eq ^ Ej. But a 
spin flip in the logical basis can also be performed exactly. 
It requires that we switch from the idle point Tyidie to the 
point where the effective magnetic field is orthogonal to 
the idle one, rj = 77idio + Tr/2. This changes the Hamil- 
tonian from H = — iAi?(ryidic)pz to H = —^AE{r]idic + 
tt/2)px. To achieve this, the dimensionless gate charge 
rig should be increased by i?j/(4i?c sin277idic)- For the 
limit discussed above, ryidio <C 1, this operating point is 
close to the degeneracy point, i] = n/2. 

(3) An alternative way to manipulate the qubit is to use 
resonant pulses, i.e., ac-pulses with frequency close to the 
qubit's level-spacing. We do not describe this technique 
here as it is well known from NMR methods. 

(4) So far we were concerned with the time dependence 
during elementary rotations. However, frequently the 
quantum state should be kept unchanged for some time, 
for instance, while other qubits are manipulated. Even 
in the idle state, 77 — jyidic, because the energies of the 
two eigenstates differ, their phases evolve relative to each 
other. This leads to the 'coherent oscillations', typical for 
a quantum system in a superposition of eigenstates. We 
have to keep track of this time dependence with high 
precision and, hence, of the time to from the very begin- 
ning of the manipulations. The time-dependent phase 
factors can be removed from the eigenstates if all the 
calculations are performed in the interaction representa- 
tion, with zero-order Hamiltonian being the one at the 
idle point. However, the price for this simplification is 
an additional time dependence in the Hamiltonian dur- 
ing operations, introduced by the transformation to the 



interaction representation. This point has been discussed 
in more detafl by Makhlin et al. (20006). 

(5) The choice of the logical basis of the qubit is by no 
means unique. As follows from the preceding discussion, 
we can perform x- and z-rotations in the charge basis, 
It), II), which provides sufficient tools for any unitary 
operation. On the other hand, since we can perform any 
unitary transformation, we can choose any other basis as 
logical basis as well. The Hamiltonian at the idle point 
is diagonal in the eigenbasis ( ^.7[ ), while the controllable 
part of the Hamiltonian, the charging energy, favors the 
charge basis. The preparation procedure (thermal relax- 
ation at the idle point) is easier described in the eigen- 
basis, while coupling to the meter (see Section ^ is di- 
agonal in the charge basis. So, the choice of the logical 
states remains a matter of convention. 

(6) A final comment concerns normal-metal single- 
electron systems. While they may serve as classical bits 
and logic devices, they are ruled out as potential quantum 
logic devices. The reason is that due to the large num- 
ber of electron states involved, their phase coherence is 
destroyed in the typical sequential tunneling processes. 



B. Charge qubit with tunable coupling 

A further step towards the ideal model (|Al|), where 
the tunneling amplitude (x-component of the field) is 
controlled as well, is the ability to tune the Josephson 
coupling. This is achieved by the design shown in Fig. ||, 
where the single Josephson junction is replaced by two 
junctions in a loop configuration (Makhlin et al. 1999). 
This dc-SQUID is biased by an external flux <i>x, which 
is coupled into the system through an inductor loop. 
If the self-inductance of the SQUID loop is low (Tin- 
kham 1996), the SQUID-controUed qubit is described by 
a Hamiltonian of the form (2.1) with modified potential 
energy: 



-£;?cos [ e + TT— 



-i?,?cos(e-7r|^) = 



- 2E'^ cos (^^^^ cos e . (2.10) 

Here $0 = hc/2e denotes the flux quanturp.. We as- 
sumed that the two junctions are identical □ with the 
same E^. The effective junction capacitance is the sum 
of individual capacitances of two junctions, in symmetric 
cases Cj = 2C?. 



*While this cannot be guaranteed with high precision in an 
experiment, we note that the effective Josephson coupling can 
be tuned to zero exactly by a design with three junctions. 
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FIG. 3. A charge qubit with tunable effective Josephson cou- 
pling. The single Josephson junction is replaced by a flux- 
threaded SQUID. The flux in turn can be controlled by a 
current carrying loop placed on top of the structure. 

When parameters are chosen such that only two charge 
states play a role, we arrive again at the Hamiltonian 
(^.3[), but the effective Josephson coupling, 



Ej{^^) = 2£;5cos TT 



(2.11) 



is tunable. Varying the external flux $x by amountSj-pf 
order $o changes the coupling between 2Ej and zero 
The SQUID-controUed qubit is, thus, described by the 



ideal single-bit Hamiltonian (Al), with field components 



B^{t) = dEci,{Vg{t)) and B^{t) = £;j($^(i)) controlled 
independently by the gate voltage and the flux. If we fix 
in the idle state Vg = Vdeg and $x = ^of^, the Hamil- 
tonian is zero, Ti.'^-^ = 0, and the state does not evolve 
in time. Hence, there is no need to control the total 
time from the beginning of the manipulations, Iq. If we 
change the voltage or the current, the modifled Hamilto- 
nian generates rotations around the z- or x-axis, which 
are elementary one-bit operations. Typical time spans 
of single-qubit logic gates are determined by the corre- 
sponding energy scales and are of order h/Ej, h/5Ech for 
X- and z-rotations, respectively. If at all times at most 
one of the fields, Bz{t) or Bx{t), are turned on, only the 
time integrals of their profiles determine the results of 
the individual operations. Hence these profiles can be 
chosen freely to optimize speed and simplicity of the ma- 
nipulations. 

The introduction of the SQUID permits not only sim- 
pler and more accurate single-bit manipulations, but it 
also allows us to control the two-bit couplings, as we will 



discuss next. Furthermore, it simplifies the measurement 
procedure, which is more accurate at i?j = (see Sec- 
tion 0). 



C. Controlled inter-qubit coupling 

In order to perform two-qubit logic gates we need to 
couple pairs of qubits together and to control the interac- 
tions between them. One possibility is to connect the su- 
perconducting boxes (i and j) directly, e.g., via a capaci- 
tor. The resulting charge-charge interaction is described 
by a Hamiltonian of the form (A2) with an Ising-type 
coupling term cx crlcl- Such a coupling allows an easy 
realization of a controUed-NOT operation. On the other 
hand, it has severe drawbacks. In order to control the 
two-bit interaction, while preserving the single-bit prop- 
erties discussed above, one needs a switch to turn the 
two-bit interaction on and off. Any externally operated 
switch, however, connects the system to the dissipative 
external circuit, thus introducing dephasing effects (see 
Section [V). They are particularly strong if the switch 
is attached directly and unscreened to the qubit, which 
would be required to control the direct capacitive inter- 
action. Therefore alternatives were explored where the 
control fields are coupled only weakly to the qubits. A 
solution (Makhlin et al. 1999) is shown in Fig. ||. All 
N qubits are connected in parallel to a common LC- 
oscillator mode which provides the necessary two-bit in- 
teractions. It turns out that the possibility to control the 
Josephson couplings by an applied fiux, simultaneously 
allows us to switch the two-bit interaction for each pair 
of qubits. This brings us close to the ideal model ( [A^ ) 
with a coupling term cx f yCr^. 
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FIG. 4. A register of many charge qubits coupled by the os- 
cillator modes in the LC-circuit formed by the inductor and 
the qubit capacitors. 



^If the SQUID inductance is not small, the fluc tuati ons of 
the flux within the SQUID renormalize the energy (E.10|). But 
still, by symmetry arguments, 
Josephson coupling vanishes. 



at <E>x — $o/2 the effective 



In order to demonstrate the mentioned properties of 
the coupling we consider the Hamiltonian of the chain 
(register of qubits) shown in Fig. ^. It is 
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£:j($x,)cose, 



$2 

21 ■ (^-l^) 



Here q denotes the total charge accumulated on the gate 
capacitors of the array of qubits. Its conjugate variable 
is the phase drop across the inductor, related to the 
flux by (f>/2'iT — $/$o- Furthermore, 



Cqb — 



(2.13) 



is the capacitance of the qubit in the external circuit. 

Depending on th e relations among the parameters the 
Hamiltonian ( ^.12 ) can be reduced. We first consider 
the situation where the frequency of the (g, $) oscillator, 
Lij^'^J = 1/ ^ NCqhL, is higher than typical frequencies of 
the qubit 's dynamics: 

AN) 



ch 



(2.14) 



In this case the oscillator modes are not excited, but 
still their virtual excitation produces an effective cou- 
pling between the qubits. To demonstrate this we elim- 
inate the variables q and $ and derive an effective de- 
scription in terms of the qubits' variables only. As a 
first step we perform a canonical transformation, q = 
q~iC\,,/Cj) E 2en, and 6, = 6, + 2^(Cqb/Cj) ($/$o), 
while $ and Ui are unchanged. This step leads to the 
new Hamiltonian (we omit the tildes): 
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We assume that the fluctuations of $ are weak 

Cqb 



(2.15) 



(2.16) 



since, otherwise, the Josephson tunneling terms in the 
Hamiltonian ( ^.15 ) are washed out (Shnirman et al. 
1997). Assuming ( 2.16| ) to be satisfied, we expand the 
Josephson terms in Eq. ( ^.15 ) up to linear terms in $. 
Then we can trace over the variables q and As a result 
we obtain an effective Hamiltonian, consisting of a sum 
of N one-bit Hamiltonians (J^^ and the coupling terms 



In spin-i notation this becomes i 

£^j(«'x»)£^j(«'xj) , ^ 
Hcoup = - 2^ — ^ ^r-^ — '^y'^y + const , (2.18) 
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where we introduced the scale 



Er = 



Cqb 



(2.19) 



The coupling Hamiltonian ( p.l8|) can be understood as 
the magnetic free energy of the current-biased inductor 
—LI^/2. This current is the sum of the contributions 
from the qubits with non-zero Josephson coupling, / oc 

E.i?}(<i>x.)sine.(xE.i?}('i>x.)'^;- 

Note that the strength of the interaction does not de- 
pend directly on the number of qubits N in the system. 
However, the frequency of the {q, <I>)-oscillator uj'j^ scales 
as The requirement that this frequency should 

not drop below typical eigenenergies of the qubit, ulti- 
mately, limits the number of qubits which can be coupled 
by a single inductor. 

The system with flux-controlled Josephson couplings 
E,]{^^i) and the interaction of the form ( 2.18 ) allows us 
to perform all necessary gate operations in a straightfor- 
ward way. In the idle state all Josephson couplings are 
turned off and the interaction (2.18) is zero. Depending 
on the choice of the idle state we may also tune the qubits 
by their gate voltages to the degeneracy points, which 
makes the Hamiltonian vanish, Ti. — Q. The interaction 
Hamiltonian remains zero during one-bit operations, as 
long as we perform one such operation at a time, i.e., 
only for one qubit we have E^ = Ej{^^i) ^ 0. To per- 
form a two-bit operation for any pair of qubits, say, i and 
J, £'} and E3^ are switched on simultaneously, yielding 
the Hamiltonian 



El 



(2.20) 



While (2.2C) is not identical to the form ( |A2| ) it equally 
well allows the relevant non-trivial two-bit operations, 
which, combined with the one-bit operations discussed 
above, provide a universal set of gates. 

A few comments should be added: 
(1) We note that typical time spans of two-bit opera- 
tions are of the order YiEl IE'\ . It follows from the con- 
ditions (2.14) and ( p. 16 ) that the interaction energy is 
never much larger than E^. Hence, at best the two-bit 
gate can be as fast as a single-bit operation. 



coup 
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(2.17) 



''While the expression (2.18) is valid only in leading order in 
an expansion in E}/hujj^Q, higher terms also vanish when the 



Josephson couplings are put to zero, 
in the idle periods persists. 



Hence, the decoupling 
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(2) It may be difficult to fabricate a nano-meter-scale 
inductor with the required inductance L, in particular, 
since it is not supposed to introduce stray capacitances. 
However, it is possible to realize such an element by a 
Josephson junction in the classical regime (with negligi- 
ble charging energy) or an array of junctions. 

(3) The design presented above does not permit perform- 
ing singe- or two-bit operations simultaneously on differ- 
ent qubits. However, this becomes possible in more com- 
plicated designs where parts of the many-qubit register 
are separated, e.g., by switchable SQUIDs. 

(4) In the derivation of the qubit interaction presented 
here we assumed a dissipation-less high-frequency oscil- 
lator mode. To minimize dissipation effects, the circuit, 
including the inductor, should be made of superconduct- 
ing material. Still at finite frequencies some dissipation 
will arise. To estimate its influence, the effect of an 
Ohmic resistance R in the circuit has been analyzed by 
Shnirman et al. (1997), with the result that the inter- 
qubit coupling persists if the oscillator is underdamped, 
R ^ ^JL/NCc^h■ In addition the dissipation causes de- 
phasing. An estimate of the resulting dephasing time 
can be obtained along the lines of the discussion in Sec- 
tion For a reasonably low-loss circuit the dephasing 
due to the coupling circuit is weaker than the influence 
of the external control circuit. 




'g2 



c 



FIG. 5. A register of charge qubits coupled to an inductor 
via separate capacitors Cl ~ Cj, independent from the gate 
capacitors Cg. 



(5) The interaction energy ( 2.18| ) involves via El the ra- 
tio of Cj and Cqb . The latter effectively screens the qubit 
from electromagnetic fluctuations in the voltage source's 
circuit, and hence should be taken as low as possible (see 
Section IV). Consequently, to achieve a reasonably high 
interaction strength and hence speed of two-bit opera- 
tions a large inductance is needed. For typical values 
of i?j ~ 100 mK and Cg/Cj ^ 0.1 one needs an in- 
ductance of L > 1 /iH in order not to have the two-bit 
operation more than ten times slower than the single-bit 
operation. However, large values of the inductance are 
difficult to reach without introducing large stray capaci- 



tances. To overcome this problem Makhlin et al. (2000a) 
suggested to use separate gate capacitors to couple the 
qubits to the inductor, as shown in Fig. |[ As long as 
the superconducting circuit of the inductor is at most 
weakly dissipative, there is no need to screen the qubit 
from the electromagnetic fluctuations in this circuit, and 
one can choose Cl as large as Cj (still larger Cl would 
decrease the charging energy Ec of the superconduct- 
ing box), which makes the relevant capacitance ratio in 
Eq. (2.17) of order one. Hence, a fairly low inductance in- 
duces an interaction of sufflcient strength. For instance, 
for the circuit parameters mentioned above, L ^ 10 nH 
would suffice. At the same time, potentially dephasing 
voltage fluctuations are screened by Cg <C Cj. 
(6) So far we discussed manipulations on time scales 
much slower than the eigenfrequency of the LC-circuit, 
which leave the LC-oscillator permanently in the ground 
state. Another possibility is to use the oscillator as a bus 
mode, similar to the techniques used for ion traps. In this 
case an ac-voltage with properly chosen frequency is ap- 
plied to a qubit to entangle its quantum state with that 
of the LC-circuit (for instance, by exciting the oscillator 
conditionally on the qubit 's state). Then by addressing 
another qubit one can absorb the oscillator quantum, 
simultaneously exciting the second qubit. As a result, 
a two-qubit unitary operation is performed. This cou- 
pling via real excitations is a first-order process, as op- 
posed to the second-order interaction (2.18). Hence, this 



method allows for faster two-qubit operations. Apart 
from this technical advantage, the creation of entangle- 
ment between a qubit and an oscillator would by itself 
be a very interesting experimental achievement (Buisson 
& Hekking 2000). 



D. Experiments with Josephson charge qubits 

Several of the concepts and properties described above 
have been verified in experiments. This includes the 
demonstration of superpositions of charge states, the 
spectroscopic verification of the spectrum, and even the 
demonstration of coherent oscillations. 

In a superconducting charge box the coherent tunnel- 
ing of Cooper pairs produces eigenstates which are gate- 
voltage dependent superpositions of charge states. This 
property has been first observed, in a somewhat indi- 
rect way, in the dissipative current through supercon- 
ducting single-electron transistors. In this system single- 
electron tunneling processes (typically treated in per- 
turbation theory) lead to transitions between the eigen- 
states. Since the eigenstates are not pure charge states, 
the Cooper-pair charge may also change in a transition. 
In the resulting combination of coherent Cooper-pair tun- 
neling and stochastic single-electron tunneling the charge 
transferred is not simply e and the work done by the 
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voltage source not simply eV. (In an expansion in the 
Josephson coupling to n-th order the charge (2n + l)e 
is transferred.) As a result a dissipative current can be 
transferred at subgap voltages. The theoretical analysis 
predicted a richly structured I-V characteristic at sub- 
gap voltages (Averin & Aleshkin 1989, Maassen van den 
Brink et al. 1991, Siewert & Schon 1996), which has been 
qualitatively confirmed by experiments (Maassen van den 
Brink et al. 1991, Tuominen et al. 1992, Hadley et al. 
1998). 

A more direct demonstration of eigenstates which arise 
as superpositions of charge states was found in the Saclay 
experiments (Bouchiat 1997, Bouchiat et al. 1998). In 
their setup (see Fig. ||) a single-electron transistor is cou- 
pled to a superconducting charge box (in the same way as 
the measurement setup to be discussed in Section 0) and 
the expectation value of the charge of the box was mea- 
sured. When the gate voltage is increased adiabatically 
this expectation value increases in a series of rounded 
steps near half-integer values of rig. At low temperatures 
the width of this transition agrees quantitati vely wit h the 
predicted ground state properties of Eqs. (2.3 2.7). At 
higher temperatures, the excited state contributes, again 
as expected from theory. 




FIG. 6. Scanning electron micrograph of a Cooper-pair box 
coupled to a single-electron transistor used in the experiments 
of the Saclay group (Bouchiat 1997, Bouchiat et al. 1998). 

Next we mention the experiments of Nakamura et al. 
(1997) who studied the superconducting charge box by 
spectroscopic means. When exposing the system to radi- 
ation they found resonances (in the tunneling current in 
a suitable setup) at frequencies corresponding to the dif- 
ference in the energy between excited and ground state, 
again in quantitative agreement with the predictions of 



The most spectacular demonstration so far of the con- 
cepts of Josephson qubits has been provided by Naka- 
mura et al. (1999). Their setup is shown in Fig. |^. In 
the experiments the Josephson charge qubit is prepared 
far from the degeneracy point for sufficiently long time 
to relax to the ground state. In this regime the ground 
state is close to a charge state, say, |t)- Then the gate 
voltage is suddenly switched to a different value. Let us, 
first, discuss the case where it is switched precisely to the 
degeneracy point. Then the initial state, a pure charge 
state, is an equal-amplitude superposition of the ground 
state |0) and the excited state |1). These two eigenstates 
have different energies, hence in time they acquire differ- 
ent phase factors: 



\^P{t)) = e-*^° |0) + e-'^i */'' |1) 



(2.21) 



After a delay time At the gate voltage is switched back 
to the original gate voltage. Depending on the delay, the 
system then ends up either in the ground state |t) [for 
{El — Eo)At/h = 2vTT with v integer], in the excited state 
\[) [for {El — Eo)At/h = {2v + l)7r], or in general in a 
At-dependent superposition. The probability that, as a 
result of this manipulation, the qubit is in the excited 
state is measured by monitoring the current through a 
probe junction. In the experiments this current was av- 
eraged over many repeated cycles, involving relaxation 
and switching processes, and the oscillatory dependence 
on At described above has been observed. 




Eq. (2.3) 



FIG. 7. Micrograph of a Cooper-pair box with flux-controlled 
Josephson junction and a probe junction (Nakamura et al. 
1999). 

In fact even more details of the theory have been quan- 
titatively confirmed. For instance one expects and finds 
an oscillatory behavior also when the gate voltage is 
switched to a point different from the degeneracy point, 
with the frequency of oscillations being a function of this 
gate voltage. Secondly, the frequency of the coherent 
oscillations depends on the Josephson coupling energy. 
The latter can be varied, since the Josephson couphng is 
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controlled by a flux-threaded SQUID (see Fig. ||). Also 
this aspect has been verified quantitatively. 

The coherent oscillations with a period of roughly 
100 psec could be observed in the experiments of Naka- 
mura et al. (1999) for at least 2 nsec 0. This puts a lower 
limit on the phase coherence time and, in fact, repre- 
sents its first direct measurement in the time domain. Es- 
timates show that a major contribution to the dephasing 
is due to the measurement process by the probe junction 
itself. In the experiments so far the detector was perma- 
nently coupled to the qubit and observed it continuously. 
Still, information about the quantum dynamics could be 
obtained since the coupling strength was optimized: it 
was weak enough not to destroy the quantum time evo- 
lution too fast and strong enough to produce a sufficient 
signal. A detector which does not induce dephasing dur- 
ing manipulations should significantly improve the op- 
eration of the device. In Section we suggest to use 
a single-electron transistor, which performs a quantum 
measurement only when switched to a dissipative state. 

So far only experiments with single qubits have been 
demonstrated. Obviously the next step is to couple two 
qubits and to create and detect entangled states. Ex- 
periments in this direction have not been successful yet, 
partially because of difficulties as, for instance, dephas- 
ing due to fiuctuating background charges. However, the 
experience with experiments with single qubits demon- 
strates that extensions to coupled qubits should be pos- 
sible as well. 

E. Adiabatic charge manipulations 

Another qubit design, based on charge degrees of free- 
dom in Josephson junction systems was proposed by 
Averin (1998). It also allows controlling the two-bit cou- 
pling at the price of representing each qubit by a chain of 
Josephson coupled islands. The basic setup is shown in 
Fig. H. Each superconducting island (with index i) is bi- 
ased via its own gate capacitor by a gate voltage Vi . The 
control of these voltages allows moving the charges along 
the chain similar to the adiabatic pumping of charges 
in junction arrays (see e.g. Pekola et al. 1999). The 
capacitances of the Josephson junctions as well as the 
gate capacitances are small enough so that the typical 
charging energy prevails over the Josephson coupling. In 
this regime the appropriate basis is that of charge states 
|ni, 712, ■ • ■), where rii is the number of extra Cooper pairs 
on island i. There exist gate voltage configurations such 
that the two charge states with the lowest energy are al- 
most degenerate, while all other charge states have much 



In later experiments the same group reported phase coher- 
ence times as long as 5 nsec (Nakamura et al. 2000) . 



higher energy. For instance, if all voltages are equal ex- 
cept for the voltages Vm and Vi at two sites, m and I, one 
can achieve the situation where the states |0, 0, 0, . . .) and 
|0, . . . , — Im, 0, . . . , 1/, . . .) are degenerate. The subspace 
of these two charge states is used as the logical Hilbert 
space of the qubit. They are coupled via the Josephson 
tunneling across the jm — ^| — 1 intermediate junctions. 




FIG. 8. Two coupled qubits as proposed by Averin (1998). 

The parameters of the qubit 's Hamiltonian can be 
tuned via the bias voltages. Obviously, the bias energy 
Bz{Vi,V2, ■■■) between these two states can be changed 
via the local voltages, Vi and Vm- Furthermore, the effec- 
tive tunneling amplitude BxiVi, V2, ■■■) can be tuned by 
adiabatic pumping of charges along the chain, changing 
their distance |to — Z| and, hence, the effective Josephson 
coupling, which depends exponentially on this distance. 
(The Cooper pair must tunnel via \m — l\ — 1 virtual 
charge states with much higher energy.) 

An inter-qubit interaction can be produced by plac- 
ing a capacitor between the edges (outer islands) of two 
qubits. If at least one of the charges in each qubit is 
shifted close towards this capacitor, the Coulomb inter- 
action leads to an interaction of the type J^^cr^cr^. The 
resulting two-bit Hamiltonian is of the form 

(2.22) 

For controlled manipulations of the qubit the coefficients 
of the Hamiltonian are modified by adiabatic motion of 
the charges along the junction array. The adiabaticity is 
required to suppress transitions between different eigen- 
states of the qubit system. While conceptually satisfying, 
this proposal appears difficult to implement: It requires 
many gate voltages for each qubit. Due to the complexity 
a high accuracy of the operations is required. Its larger 
size as compared to simpler designs makes the system 
more vulnerable to dephasing effects, e.g., due to fiuctu- 
ations of the offset charges. 

Adiabatic manipulations of the Josephson charge qubit 
can lead to Berry phases. Falci et al. (2000) suggested 
that a Berry phase can accumulate during suitable ma- 
nipulations of a fiux controlled charge qubit with an 
asymmetric dc-SQUID, and that it can be detected in 
an experiment similar to that of Nakamura et al. (1999). 
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If the bare Josephson couplings of the SQUID loop are 
Ej and the effective Josephson energy is given by (cf. 
Eq.(|lOl)j 



-E] cos ( e + TT— ) - £:? cos I e - TT — 



(2.23) 



Hence, the corresponding Hamiltonian of the qubit has 
all three components of the effective magnetic field: — 
(S] + S2)cos(7r$x/$o) andSy = (S^-S]) sin(7r$,/$o), 
while Bz is given by Eq.(|]J). With three non-zero field 
components, adiabatic changes of the control parameters 
Vg and <i>x may result in B enclosing a non-zero solid 
angle. This results in a Berry phase shift, 7b, between 
the ground and excited states. In general, a dynamic 
phase / AE{t)dt is also accumulated in the process. To 
single out the Berry phase, Falci et al. (2000) suggested 
to encircle the loop in parameter space back and forth, 
with a NOT operation performed in between. The latter 
exchanges the ground and excited state, and, as a result, 
the dynamic phases accumulated during both paths can- 
cel. At the same time the Berry phases add up to 273. 
This phase shift can be measured by a procedure similar 
to that used by Nakamura et al. (1999): the system is pre- 
pared in a charge state away from degeneracy, abruptly 
switched to the degeneracy point where adiabatic ma- 
nipulations and the NOT gate are performed, and then 
switched back. Finally, the average charge is measured. 
The probability to find the qubit in the excited charge 
state, sin^ 27B, reflects the Berry phase. 

The experimental demonstration of topological phases 
in Josephson junction devices would constitute a new 
class of "macroscopic" quantum effects in these systems. 
They can be performed with a single Josephson qubit 
in a design as used by Nakamura et al. (1999) and thus 
appear feasible in the near future. 



III. QUBITS BASED ON THE FLUX DEGREE OF FREEDOM 

In the previous Section we described the quantum dy- 
namics of low capacitance Josephson devices where the 
charging energy dominates over the Josephson energy, 
Ec ^ Ej , and the relevant quantum degree of freedom 
is the charge on superconducting islands. We will now re- 
view the quantum properties of superconducting devices 
in the opposite regime, Ej ^ Eq, where the flux is the 
appropriate quantum degree of freedom. These systems 
were proposed by Caldeira & Leggett (1983) in the mid 
80s as test objects to study various quantum mechanical 
effects. This includes the 'macroscopic quantum tunnel- 
ing' of the phase (or flux) as well as resonance tunneling. 
Both had been observed in several experiments (Voss & 
Webb 1981, Martinis et al. 1987, Clarke et al. 1988, Rouse 
et al. 1995). Another important quantum effect has been 
reported recently: The groups in Stony Brook (Fried- 
man et al. 2000) and in Delft (van der Wal et al. 2000) 



demonstrated in experiments the avoided level crossing 
due to coherent tunneling of the ffux in a double well po- 
tential. In principle, all other manipulations discussed in 
the previous section should be possible with Josephson 
flux devices as well. They have the added advantage not 
to be sensitive to fluctuations in the background charges. 
However, attempts to observe 'macroscopic quantum co- 
herent oscillations' in Josephson flux devices have not 
been successful yet (Leggett 1987, Tesche 1990). 



A. Josephson flux (persistent current) qubits 

We consider superconducting ring geometries inter- 
rupted by one or several Josephson junctions. In these 
systems persistent currents flow and magnetic fluxes are 
enclosed. The simplest design of these devices is an rf- 
SQUID, which is formed by a loop with one junction, as 
shown in Fig. || a. The phase difference across the junc- 
tion is related to the flux <I> in the loop (in units of the 
flux quantum $0 = h/2e) by (p/2Tr — $/$o + integer. An 
externally applied flux <I>x biases the system. Its Hamil- 
tonian, with Josephson coupling, charging energy, and 
magnetic contributions taken into account, thus reads 



n 



-Ej cos 2tt 



$0 



($--$x) 

2L 



2Cj 



(3.1) 



Here L is the self-inductance of the loop and Cj the ca- 
pacitance of the junction. The charge Q = —ihd/d^ on 
the leads is canonically conjugate to the flux <I>. 




(J) 



X 



FIG. 9. a) The rf-SQUID, a simple loop with a Josephson 
junction, forms the simplest Josephson flux qubit. b) An im- 
proved design of a flux qubit. The flux $x in the smaller loop 
controls the effective Josephson coupling of the rf-SQUID. 

If the self-inductance is large, such that = 
i?j/($Q/47r^L) is larger than 1 and the externally ap- 
plied flux <I>x is close to $o/2, the first two terms in the 
Hamiltonian form a double- well potential near <f> = $o/2. 
At low temperatures only the lowest states in the two 
wells contribute. Hence the reduced Hamiltoni an o f 
this effective two-state system again has the form (^.5|), 
Wctri = —\Bz^z ~ \Bx(yx- The diagonal term B^ is the 
bias, i.e., the asymmetry of the double well potential, 
given for /3l — 1 <C 1 by 



S.($x) = 4^v/6(/3l - 1) Ei ($x/$o - 1/2) . (3.2) 
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It can be tuned by the applied flux $x- The off-diagonal 
term describes the tunneling amplitude between the 
wells. It depends on the height of the barrier and thus 
on Ej. This Josephson energy, in turn, can be controlled 
if the junction is replaced by a dc-SQUID, as shown in 
Fig. P b, introducing the flux as another control vari- 
able B. With these two external control parameters the 
elementary single-bit operations, i.e., z- and x-rotations 
can be performed, equivalent to the manipulations de- 
scribed for charge qubits in the previous section. Also 
for flux qubits we can either perform the operations by 
sudden switching of the external fluxes $x and $x for a fi- 
nite time, or we can use ac-fields and resonant pulses. To 
enable coherent manipulations the parameter should 
be chosen larger than unity (so that two wells with well- 
defined levels appear) but not much larger, since the re- 
sulting large separation of the wells would suppress the 
tunneling. 

The rf-SQUID described above had been discussed in 
the mid 80s as a realization of a 2-state quantum system. 
Some features of macroscopic quantum behavior were 
demonstrated such as 'macroscopic quantum tunneling' 
(MQT) of the flux, resonant tunneling and level quanti- 
zation (Voss & Webb 1981, Martinis et al. 1987, Clarke 
et al. 1988, Rouse et al. 1995, Silvestrini et al. 1997). 
However, only very recently coherent superpositions of 
macroscopically different flux states have been demon- 
strated (Friedman et al. 2000, van der Wal et al. 2000). 

The group in Stony Brook (Friedman et al. 2000) 
probed spectroscopically the superposition of excited 
states in different wells. The rf-SQUID used had self- 
inductance L = 240 pH and = 2.33. A substantial 
separation of the minima of the double-well potential (of 
order <I>o) and a high inter- well barrier made the tunnel 
coupling between the lowest states in the wells negligi- 
ble. However, both wells contain a set of higher local- 
ized levels - under suitable conditions one state in each 
well - with relative energies also controlled by $x and 
$x- Being closer to the top of the barrier these states 
mix more strongly and form eigenstates, which arc su- 
perpositions of localized flux states from different wells. 
External microwave radiation was used to pump the sys- 
tem from a well-localized ground state in one well to one 
of these eigenstates. The energy spectrum of these levels 
was studied for different biases <I>x, $x, and the prop- 
erties of the model (3.1) were confirmed. In particular, 
the level splitting at the degeneracy point indicates a su- 
perposition of distinct quantum states. They differ in a 
macroscopic way: the authors estimated that two super- 
imposed flux states differ in flux by $o/4, in current by 
2-3 /iA, and in magnetic moment by 10^*^ /ie- 



®See (Mooij et al. 1999) for suggestions how to control <l>x 
independent of $x. 



The Delft group (van der Wal et al. 2000) performed 
microwave spectroscopy experiments on a similar but 
much smaller three-junction system described below. 
The ground states in two wells of the Josephson poten- 
tial landscape were probed. The obtained results verify 
the spectrum of the qubit and the level repulsion at the 
degeneracy point expected from the model Hamiltonian 
( p. 3D with the parameters B^, Bz calculated from the 
potential ( |3.3| ). Similar to the experiments of Friedman 
et al. (2000), this provides clear evidence for superposi- 
tions of macroscopically distinct phase states. 

In spite of this progress, attempts to observe 'macro- 
scopic quantum coherence', i.e., the coherent oscilla- 
tions of a quantum system prepared in a superposition 
of eigenstates have not been successful so far (Leggett 
1987, Tesche 1990). A possible reason for this failure has 
been suggested recently by Mooij et al. (1999). They 
argue that for the rf-SQUID designs considered so far 
the existence of the double-well potential requires that 
/3l > 1 which translates into a sufficiently high prod- 
uct of the critical current of the junction and its self- 
inductance. In practice, only a narrow range of circuit 
parameters is useful, since high critical currents require 
a relatively large junction area resulting in a high capaci- 
tance which suppresses tunneling. A high self-inductance 
of the rf-SQUID can be achieved only in large loops. This 
makes the system very susceptible to external noise. 

To overcome this difficulty Mooij et al. (1999) and 
Feigelman et al. (2000) proposed to use a smaller su- 
perconducting loop with three or four junctions, respec- 
tively. Here we discuss the 3-junction circuit shown 
in Fig. |l^ a,c. In this low-inductance circuit the flux 
through the loop remains close to the externally applied 
value, <I> = <I>x. Hence the phase differences across the 
junctions are constrained by fi+(p2 + f3 — 27r<I'x/*I'0j 
leaving tpi and ip2 as independent dynamical variables. 
In the plane spanned by these two variables the Joseph- 
son couplings produce a potential landscape given by 

U{ipi,ip2) = -Ej cos ipi - Ej cos ip2 

-Ej cos(27r$x/*o -m- ^2) . (3.3) 

If Ej/Ej > 0.5, a double- well potential is formed within 
each 27r x 27r cell in the phase plane. For an optimal value 
of Ej/E;i « 0.7-0.8 the cells are separated by high barri- 
ers, while tunneling between two minima within one cell 
is still possible. The lowest states in the wells form a two- 
state quantum system, with two different current config- 
urations. Mooij et al. (1999) and Orlando et al. (1999) 
discuss junctions with Ej 2 K and Ej /Eq ^ 80 and 
loops of micrometer size with very small self-inductance 
1/ ~ 5 pH (which can be neglected when calculating the 
energy levels). Typical qubit operation parameters are 
the level splitting Bz ~ 0.5 K and the tunneling am- 
plitude Bx ^ 50 niK. For the optimal choice of Ej/Ej 
the two minima differ in phases by an amount of order 
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7r/2. Due to the very low inductance and the relatively 
low critical current Ic ^ 200 nA this translates into a 
flux difference of 5$ ~ Lie 10^'^$o- While this corre- 
sponds to a still 'macroscopic' magnetic moment of 10^ to 
10^ Hb, the two basis states are similar enough to make 
the coupling to external fluctuating flelds and hence the 
dephasing effects weak (for a further discussion, see Sec- 
tion . In this respect the new design is qualitatively 
superior to the simple rf-SQUID. 



Ej 



4) 



E,«l.) 




action energy is of order MI^ where M is the mutual 
inductance and Ic = (27r/$o)^'J is the critical current 
in the junctions. For their design, Mooij et al. (1999) 
estimate the typical interaction energy to be of order 
0.01£^j ~ 50 mK in frequency units, i.e., of the order 
of single-qubit energies. For a typical rf-SQUID (Fried- 
man et al. 2000) this coupling can be even stronger than 
the tunneling rate between the flux states of the SQUID. 

In the simplest form this interaction is always turned 
on. To turn it off completely, one needs a switch which 
has to be controlled by high-frequency pulses. The re- 
lated coupling to the external circuit leads to decoher- 
ence (see the discussion at the end of this Section). An 
alternative is to keep the interaction turned on constantly 
and use ac driving pulses to induce coherent transitions 
between the levels of the two-qubit system (cf. Shnir- 
man et al. 1997, Mooij et al. 1999). A disadvantage of 
this approach is that permanent couplings result in an 
unwanted accumulation of relative phases between the 2- 
qubit states even in the idle periods. Keeping track of 
these phases, or their suppression by repeated refocus- 
ing pulses (see Section iV), require a high precision and 
complicate the operation. 
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FIG. 10. a) and c) A 3-junction loop as a flux qubit (Mooij 
et al. 1999). The reduced size and lower inductance of this sys- 
tem as compared to earlier designs (e.g., Fig. ^ a) reduce the 
coupling to the external world and hence dephasing effects, 
b) Multi-junction flux qubit with a controlled Josephson cou- 
pling (Mooij et al. 1999). Control over two magnetic fluxes, $ 
and allows one to perform all single-qubit logic operations. 



B. Coupling of flux qubits 

In order to couple different flux qubits one can use 
a direct inductive coupling (Mooij et al. 1999, Orlando 
et al. 1999) as shown by the dashed line in Fig. |l^. A 
mutual inductance between the qubits can be established 
in different ways. The dashed loop shown in the flgure 
couples the currents and fluxes in the lower parts of the 
qubits. Since fluxes through these loops control the bar- 
rier heights of the double- well potentials, this gives rise 
to the interaction term cx o'jJ.CT^. Placing the loop differ- 
ently produces in addition contributions to the interac- 
tion Hamiltonian of the form o'^o'^. The typical inter- 
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FIG. 11. Flux qubits coupled in two ways. The dashed line 
induces a direct inductive coupling. Alternatively, an inter- 
qubit coupling is provided by the LC-circuit indicated by a 
solid line. 

A controllable inter-qubit coupling without additional 
switches is achieved in the design shown by the solid line 
in Fig. (Makhlin et al. 2000c). The couphng is me- 
diated by an LC-circuit, with self- inductance Lose and 
capacitance Cose, which is coupled inductively to each 
qubit. Similar to the design of the charge qubit register 
in Section II. C, the coupling depends on parameters of 
of individual qubits and can be controlled in this way. 
The effective coupling can be found again by integrat- 
ing out the fast oscillations in the LC-circuit. It can be 
understood in a simple way by noting that in the limit 
Cose the qubits establish a voltage drop across the 
inductor, V — M^i/L, and the Hamiltonian for the 
osciUator mode is Hose = *i>^/2Losc + Q^/2Cose - VQ, 
with the charge Q being conjugate to the flux $ through 
the LC-circuit. Here $i is the flux in the loop of qubit i, 
L is the self-inductance of the loop and M is its mutual 
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inductance with the LC-circuit. Continuing as described 



in Section 11. C we obtain the inter-qubit interaction term 
— In the hmit of weak couphng to the LC- 
circuit, we have $i = j-[Hi,^i] = S^iBlay/h, where 
is the separation between two minima of the potential 
and Bl is the tunnehng ampHtude. Hence, the interac- 
tion is given by 



Tin 



i<3 



(3.4) 



To turn off the interaction one should suppress the 
tunneling amplitudes B\.. This can be done with ex- 
ponential precision by increasing the height of the po- 
tential barrier via Note, that in this case also un- 
wanted fluctuations of B\. and resulting dephasing ef- 
fects are exponentially suppressed. All needed single and 
two-qubit manipulations can be performed by turning 
on the fields B]. and i?* , in complete analogy to what 
we discussed in Section II. C , We also encounter the 
equivalent drawbacks: the design shown in Fig. ^ does 
not allow simultaneous manipulations on different qubit 
pairs, and the conditions of high oscillator frequencies 
and weak renormalization of qubit parameters by the 
coupling, similar to Eqs.(2.14), (2.16), limit the two-qubit 
coupling energy. The optimization of this coupling re- 
quires Lose I Cose ~ RK{6^l^of{M/Lf and ljlc not 
far above the qubit frequencies. For rf-SQUIDs (Fried- 
man et al. 2000) the resulting coupling can reach the 
same order as the single-bit terms. On the other hand, 
for the design of Mooij et al. (1999), where two basis 
phase states differ only slightly in their magnetic proper- 
ties, the coupling term is much weaker than the single-bit 
energies. 



currents in the Josephson junctions while external loops 
carry no current. As a result the coupling of the qubit to 
the electromagnetic environment is substantially reduced 
and coherence is preserved longer. 

The mentioned designs are similar and we discuss them 
in parallel. loffe et al. (1999) suggested to use an SD tun- 
nel junction with the s-wave lead matched to the (110) 
boundary of the d-wave superconductor. In this geom- 
etry the first harmonic, oc cos 93, in the Josephson cou- 
pling vanishes due to symmetry reasons, and one obtains 
a bistable system with the potential energy cos 2^5 and 
minima at ±(^0 with Lp^ — 7r/2. (A similar current-phase 
relation was observed recently by Il'ichev et al. (1998) in 
a DD junction with a mismatch angle of 45°.) Zagoskin 
(1999) proposed to use DND, or D-grain-boundary-D 
Josephson junctions formed by two d-wave superconduc- 
tors with different spatial orientation of the order param- 
eter, connected by a normal metal. The energy-phase 
relation for such junctions also has two degenerate min- 
ima, at the phase differences ±c^o- The separation 2(po 
of these minima, and hence the tunneling amplitude, are 
controlled by the mismatch angle of the d-wave leads. 

In a later development a 'macroscopic analogue' of 
d-wave qubits was discussed (Blatter et al. 1999). In- 
stead of an SD-junction, it is based on a five-junction 
loop, shown in Fig. |l^, which contains one strong tt- 
j unction and four ordinary junctions. The presence of 
the TT-junction is equivalent to magnetically biasing the 
loop with a half superconducting flux quantum. Four 
other junctions, frustrated by the 7r-phase shift, have two 
lowest-energy states with the phase difference of ±7r/2 
between the external legs in the figure. In this respect the 
5-junction loop is similar to the SD-junction discussed 
above and can be called a 7r/2-junction. 



C. "Quiet" superconducting phase qubits 

The circuits considered so far in this Section are vul- 
nerable to external noise. First, they need for their op- 
eration an external bias in the vicinity of $o/2, which 
should be kept stable for the time of manipulations. In 
addition, the two basis flux states of the qubit have differ- 
ent current configurations, which may lead to magnetic 
interactions with the environment and possible crosstalk 
between qubits. To a large extent the latter effect is 
suppressed already in the design of Mooij et al. (1999). 
To further reduce these problems several designs of so- 
called "quiet" qubits have been suggested (loffe et al. 
1999, Blatter et al. 1999, Zagoskin 1999, Blais & Zagoskin 
2000). They are based on intrinsically doubly degener- 
ate systems, e.g., Josephson junctions with d-wave leads 
and energy-phase relation (e.g., cos 20) with two minima, 
or the use of 7r-j unctions which removes the need for a 
constant magnetic bias near $0/2. As a result the rele- 
vant two states differ only in their distribution of internal 
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FIG. 12. A five-junction loop, a basic bistable element of a 
"quiet" superconducting qubits (Blatter et al. 1999), is made 
of four ordinary junctions and one stronger 7r-junction. In two 
stable configurations the phase difference across this element 
is ±7r/2. 

In all these designs the bistability is a consequence of 
the time-reversal symmetry (which changes the signs of 
all the phases) of the Hamiltonian. Thus the degener- 
acy persists also in systems containing different Joseph- 
son junctions, although the phase differences in the two 
lowest-energy states and their separation can change. If 
charging effects with Eq <C i?j are taken included, one 
arrives at a double-well system with tunneling between 
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the wells. Such a qubit can be operated by connecting 
or disconnecting it from external elements, as described 
below. 

The first issue to be addressed is the question how to 
store the qubit's state, i.e., how to freeze the evolution. 
This can be achieved by connecting the qubit in parallel 
to a large capacitor (lofFe ct al. 1999). This makes the 
phase degree of freedom very massive, thus suppressing 
the tunneling and restoring the needed degeneracy. In 
order to perform a (Tj;-rotation the inter-well tunneling is 
turned on by disconnecting the capacitor. This means a 
switch is needed in the circuit. 

The (T^-rotation or phase shift can be accomplished by 
lifting the degeneracy between the wells. This can be 
done by connecting another, much stronger '7r/2-junction 
(DS-junction or 5-junction loop) and a weak ordinary, s- 
wave junction (with Josephson energy (x cos ip) in series 
to the qubit, to form a closed loop. This again requires a 
switch. The auxiliary 7r/2-junction shifts the phase dif- 
ferences of the potential minima of the qubit to and 
TT. Hence the s-junction is in the ground state or frus- 
trated depending on the qubit's phase drop. The cor- 
responding energy difference produces the needed phase 
shift between two qubit's states. 

To perform two-qubit manipulations and control the 
entanglement loffe et al. (1999) proposed to form a loop, 
connecting in series two qubits and one s-junction with 
weak Josephson coupling E"^ <^ E_t. The phase state 
of each qubit is characterized by the phase difference of 
±ipo, i-e., the total phase drop on the qubits is equal to 
±2(^0 or depending on whether the qubits are in the 
same state or in a different ones. When the connection 
between the qubits and the s-junction is turned on, this 
phase drops across the s-junction. and its energy differs 
by -Ej(l — cos2(^o) for the states |00), |11) as compared 
to the states |01), |10). The net effect is an Ising-type 
interaction between the pair of qubits, which allows per- 
forming unitary two-qubit transformations. 

Another way of operation was discussed by Blais & 
Zagoskin (2000). They suggested to use a magnetic force 
microscope tip for single-bit manipulations (local mag- 
netic field lifts the degeneracy of two phase states) and 
for the read-out of the phase state. The tip should be 
moved towards or away from the qubit during manipula- 
tions. The short time scales of qubit operation make this 
proposal difficult to reahze. 

Even in "quiet" designs, in both SD and DD systems, 
there are microscopic persistent currents flowing inside 
the junctions which differ for the two logical states (Blat- 
ter et al. 1999, Zagoskin 1999). These weak currents still 
couple to the outside world and to other qubits, thus 
spoiling the ideal behavior. Furthermore, all the designs 
mentioned require externally operated switches to con- 
nect and disconnect qubits. We discuss the associated 
problems in the following subsection. 



To summarize, the "quiet" designs require rather com- 
plicated manipulations as well as circuits with many junc- 
tions, including 7r-junctions or d-wave junctions, which 
are difficult to fabricate in a controlled and reliable way. 
In addition, many constraints imposed on the circuit pa- 
rameters (in particular, on the hierarchy of Josephson 
couplings) appear difficult to satisfy. In our opinion the 
"quiet" phase qubit designs belong to a higher complex- 
ity class than the previously discussed charge and flux 
qubits, and their experimental realization may remain a 
challenge for some time. 

D. Switches 

Switches may be used in variety of contexts in quantum 
nano-circuits. They are needed, e.g., for a direct capaci- 
tive coupling between charge qubits or magnetic coupling 
of flux qubits. They are also a major tool for controlling 
the dynamics of "quiet" qubits. Ideal switches should 
decouple qubits from the environment and at the same 
time let through control signals. They should operate on 
the very fast time scale of the qubit dynamics and have a 
high "switching ratio" , that is the ratio of the interaction 
with the switch in the on- or off-state. Such switches are 
hard to realize. In this subsection we compare the char- 
acteristics of several Josephson-junction-based switches 
and the associated problems. 

Possible switches are dc-SQUIDs as well as SSETs 
(single-Coopcr-pair transistors) in a mode where they act 
as Josephson junctions with an externally controlled cou- 
pling. Then the switching ratio is the ratio of the minimal 
and maximal values of the coupling. In a dc- SQUID with 
Josephson energies of its junctions equal to Ej and Ej, 
this ratio is (E] ~ Ej)/{E] + Ej). It reached a value 
below 1% in the experiment of Rouse et al. (1995). How- 
ever, fast switching of the bias flux may be difficult to 
perform. In a SSET the effective coupling is controlled 
by a gate voltage, which can be switched fast. However, 
the switching ratio of order Ej / Eq {Ej and Ec are char- 
acteristics of the SSET) is hardly below several percent. 
These non-idealities lead to unwanted interactions when 
the switch is supposed to be disconnected. 

Since a dc-SQUID requires an external bias to be oper- 
ated as a switch, Blatter et al. (1999) suggested a similar 
construction with the bias provided by 7r/2-junctions in- 
stead of an external magnetic field. Namely, one can 
insert two tt/ 2- junctions into one arm of the SQUID- 
loop. Depending on whether the phase drops across 
these junctions are equal or opposite, they simulate an 
external bias of a half flux quantum or no bias. Accord- 
ingly, the Josephson couplings of two s-j unctions in the 
SQUID add up or cancel each other. The switching is re- 
alized via a voltage pulse which drives one 7r/2-junction 
between +tt/2- and — 7r/2-states. Blatter et al. (1999) 
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also suggest to use an array of n such switches, reducing 
the overall Josephson coupling in the off-state by a fac- 
tor [(£'] — E'j)/Ec\"-. Unfortunately, in the on-state the 
overall coupling through the array is also reduced with 
growing n, although this reduction may be weaker than 
in the off-state, i.e., the switching ratio increases with n. 
Still the quality of the switch in the on-state is reduced. 
Moreover, to operate the switch one would need to send 
simultaneously voltage pulses to n intermediate elements 
which complicates the operation. Note that this design 
is reminiscent of the qubit design proposed by Averin 



(1998) which is presented in Subsection II. E. However 



while Blatter et al. (1999) suggest to control the cou- 
pling, cx (E'j/E'c)", by controlling Averin proposes 
to change the distance n of the tunneling process. 

While switches of the type as described above may be 
useful in first experiments with simple quantum nano- 
circuits, further work is needed before they can be used 
in more advanced designs which require high precision of 
manipulations and phase coherence for a long time. 



IV. ENVIRONMENT AND DISSIPATION 
A. Identifying the problem 

For an ideal quantum system the time evolution is de- 
scribed by deterministic, reversible unitary operations. 
The concepts of quantum state engineering and compu- 
tation heavily rely on this quantum coherence, with many 
potential applications requiring a large number of coher- 
ent manipulations of a large number of qubits. On the 
other hand, for any real physical quantum system the 
time evolution may be disturbed in various ways, and 
the number of coherent manipulations is limited. Possi- 
ble sources of errors are inaccuracies in the preparation of 
the initial state, inaccuracies in the manipulations (logic 
gates), uncontrolled couplings between qubits, undesired 
excitations out of the two-state Hilbert space (Fazio et al. 
1999), and - unavoidable in devices which are to be con- 
trolled externally - interactions with the environment. 
Due to the coupling to the environment the quantum 
state of the qubits gets entangled with the environmental 
degrees of freedom. As a consequence the phase coher- 
ence is destroyed after a time scale called the dephasing 
time. In this Section we will describe the influence of 
the environment on the qubit. We determine how the 
dephasing time depends on system parameters and how 
it can be optimized. 

Some of the errors can be corrected by 'software' tools. 
One known from NMR and, in particular, NMR-based 
quantum logic operations (see e.g. Chang 1998) are the 
refocusing techniques. They serve to suppress the effects 
of undesired terms in the Hamiltonian, e.g., deviations 
of the single-bit field terms from their nominal values 
or uncontrolled interactions like stray direct capacitive 



couplings of charge qubits or inductive couplings of flux 
qubits. As an example we consider the error due to a 
single-bit term 5Bx(Jx, which after some time has pro- 
duced an unwanted rotation by a. Refocusing is based 
on the fact that a 7r-pulse about the z-axis reverses the 
influence of this term, i.e., Uz{TT)Ux{a)Uz{TT) = Ux{—Oi). 
Hence, fast repeated inversions of the bias Bz{t) (with 
\Bz\ ^ SBx) eliminate the effects associated with SBx- 
The technique can also be applied to enhance the pre- 
cision of non-ideal control switches: one first turns off 
the coupling term to a low value and then further sup- 
presses it by refocusing. The examples demonstrate that 
refocusing requires very fast repeated switchings with a 
period much shorter than the elementary operation time. 
This may make it hard to implement. 

It was therefore a major breakthrough when the con- 
cepts of quantum error-correction were discovered (see 
e.g. Steane 1998, PreskiU 1998). When apphed they 
should make it possible, even in the presence of dephas- 
ing processes - provided that the dephasing time is not 
too short - to perform coherent sequences of quantum 
manipulations of arbitrary length. The price to be paid 
is an increase in system size (by roughly an order of mag- 
nitude), and a large number of steps are needed for er- 
ror correction before another computational step can be 
performed (increasing the number of steps by roughly 3 
orders of magnitude). This imposes constraints on the 
dephasing time. The detailed analysis shows that error 
correction can be successful if the dephasing time is of 
the order of 10'* times longer than the time needed for an 
elementary logic gate. 

In the Josephson junction systems, discussed here, the 
environment is usually composed of resistive elements in 
the circuits needed for the manipulations and the mea- 
surements. They produce voltage and current noise. In 
many cases such fluctuations are Gaussian distributed 
with a Johnson-Nyquist power spectrum, coupling lin- 
early to the quantum system of interest. They can thus 
be described by a harmonic oscillator bath with suitable 
frequency spectrum and coupling strength (Leggett et al. 
1987, Weiss 1999). For charge qubits, for instance, fluc- 
tuations in the gate voltage circuit, coupling to az, as 
well as the fluctuations in the current, which control the 
Josephson energy and couple to ax, can be described 
in this way (Shnirman et al. 1997). In this section we 
will first describe these noise sources and the dephas- 
ing introduced in this way. We later comment on other 
noise sources such as telegraph noise, typically with a 
1// power spectrum due to switching two- level systems 
(e.g., 'background charge' fluctuations), or the shot noise 
resulting from the tunneling in a single-electron transis- 
tor coupled to a qubit for the purpose of a measurement. 

Depending on the relation between typical frequen- 
cies of the coherent (Hamiltonian) dynamics and the de- 
phasing rates we distinguish two regimes. In the first. 
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"Hamiltonian-dominated" regime, where the controlled 
part of the qubit Hamiltonian Hctri = — (l/2)Bc, gov- 
erning the deterministic time evolution and logic gates, 
is large, it is convenient to describe the dynamics in the 
eigenbasis of Tictri- The coupling to the environment 
is weak, hence the environment-induced transitions are 
slow. One can then distinguish two stages: a) dephasing 
processes, in which the relative phase between the eigen- 
states becomes random; h) energy relaxation processes, 
in which the occupation probabilities of the eigenstates 
change. 

In the other, "environment-dominated" regime TYctri is 
too weak to support its eigenstates as the preferred ba- 
sis. The qubit's dynamics in this situation is governed 
by dissipative terms and depends on details of the struc- 
ture of the coupling to the environment. In general the 
evolution is complicated, and the distinction between re- 
laxation and dephasing may be impossible. 

Both regimes may be encountered during manipula- 
tions. Obviously, the Hamiltonian should dominate when 
a coherent manipulation is performed. On the other 
hand, if in the idle state the Hamiltonian vanishes (a 
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{Xl) EE - ({X{t),Xit')})^ - hJiij) coth -— , (4.5) 



very useful property as outlined in Section II. A , II. B ), 
the environment-dominated regime is realized. One has 
to ensure that the phase coherence rate in this regime is 
still low enough. 



B. Spin-boson model 

Before we proceed discussing specific physical sys- 
tems, we recall what is known about the spin-boson 
model, which has been studied extensively (see reviews 
by Leggett et al. 1987, Weiss 1999). It models the envi- 
ronment as an oscillator bath coupled to one component 
of the spin. The Hamiltonian reads 



(4.1) 



where the bath spectral density is defined by 

J(c.)^5^^5(c.-..J. 
2 ^ mau;a 



(4.6) 



This spectral density has typically a power-law behavior 
at low frequencies (Leggett et al. 1987). Of particular in- 
terest is Ohmic dissipation, corresponding to a spectrum 



J{Lij) — — ahuj 



(4.7) 



which is linear at low frequencies up to some high- 
frequency cutoff uJc- The dimensionless parameter a re- 
flects the strength of dissipation. Here we concentrate on 
weak damping, a <C 1, since only this regime is relevant 
for quantum state engineering. But still the Hamiltonian- 
dominated and the environment-dominated regimes are 
both possible depending on the ratio between the energy 
scale AE = \J Bl + -B^ , characterizing the coherent evo- 
lution, and the dephasing rate (to be determined below). 

The Hamiltonian-dominated regime is realized when 
AE 3> akBT. In this regime it is natural to describe 
the evolution of the system in the eigenbasis (2.7) which 
diagonalize Tictri: 

n = -^AEp, + (sin?7 + COST? p,) X + Hb ■ (4.8) 

Two different time scales characterize the evolu- 
tion (Weiss & WoUensak 1989, Gorlich et al. 1989, Weiss 
1999). On a first, dephasing time scale the off- 
diagonal (in the preferred eigenbasis) elements of the 
qubit's reduced density matrix decay to zero. They are 
represented by the expectation values of the operators 
p± = (1/2) (pa; ± ipy). Dephasing leads to the following 
time-dependence (at long times): 



{p±{t)) ^ {p±m e 



(4.9) 



where 
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(4.2) 
(4.3) 



is the controlled part of the Hamiltonian (cf. Eqs. 2.3 
and 2.6), while 
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On the second, relaxation time scale Trciax the diagonal 
entries tend to their thermal equilibrium values: 

{pM) = P.(oo) + [pM - p,{^)] e-*/^-'- , (4.10) 

where p^(oo) = tanh(A£:/2A:BT). 

The dephasing and relaxation times were originally 
evaluated for the spin boson model in a path integral 
technique (Leggett et al. 1987, Weiss 1999). The rates 
are □ 



is the Hamiltonian of the bath. The bath operator X — 
X^a ^o.Xa couples to . lu thermal equilibrium one finds 
for the Fourier transform of the symmetrized correlation 
function of this operator 



^Note that in the literature usually the evolution of {ozit)) 
as heen s tudied. To establish the m nnertinn to the results 



rinertinri 

( [l.ll|j4.12D one has to substitute Eqs. ( ^ , [l.lO| ) into the iden 
tity (Jz = cos rj pz -\- sin ri px Furthermore, we neglect renor- 
malization effects, since they are weak for a <^1. 
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coth 
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(4.11) 
(4.12) 



In some cases these results can be derived in a simple 
way, which we present here to illustrate the origin of dif- 
ferent terms. As is apparent from the Hamiltonian (4.S) 
the problem can be mapped on the dynamics of a spin- 
1/2 particle in the external magnetic field AE pointing 
in z-direction and a fluctuating field in the x-z-plane. 
The x-component of this fluctuating field, with magni- 
tude proportional to sin 77, induces transitions between 
the eigenstates (^) of the unperturbed system. Apply- 
ing the Golden Rule for this term one obtains readily the 
relaxation rate (4.11). 

The longitudinal component of the fluctuating field, 
proportional to cos?] does not induce relaxation pro- 
cesses. It does, however, contribute to dephasing since 
it leads to random fluctuations of the eigenenergies and, 
thus, to a random relative phase between the two eigen- 
states. As an example we analyze its effect on the de- 
phasing rate in an exactly solvable limit. 

The unitary operator 



C/ ^ exp ( with $ ^ ^ |^ (4.13) 
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transforms the Hamiltonian (4.1-4.4) to a rotating spin- 
frame (Leggett ct al. 1987): 

n = UnU-^ = - (1/2) A£; cos 77 cr^ 

-(l/2)A£;sin?7 (cr+e^^* + h.c.) + Hb ■ (4.14) 

Here we recognize that in the limit 77 = the spin and 
the bath are decoupled, which allows an exact treat- 
ment. The trivial time evolution in this frame, a± {t) = 
exp{^iAEt)a±{Q), translates in the laboratory frame to 



a±(0 = e^^*'*'e±'^*(°)eT^^^V±(0) . 



To average over the bath we need the correlator 



(4.15) 



P{t) = (e'*(*) e-'*(")) (e-**(*) e**^")) (4.16) 

which was studied extensively by many authors (Leggett 
et al. 1987, Panyukov & Zaikin 1988, Odintsov 1988, 
Nazarov 1989, Devoret et al. 1990). It can be expressed 
as P{t) — exp[K{t)], where 
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coth 
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(cosLut — 1) — 7 sinujt 



(4.17) 



For the Ohmic bath (|4/7|) for t > h/ik^T one has 
ReK{t) w ~{2kBT/h)Trat. Thus we reproduce Eq. (O) 



with Tip given by ( 4.12| ) in the limit 77 = 0. While it is not 
so simple to derive the general result for arbitrary 77, it is 
clear from Eqs. (4.11) and (4.12) that the effects of the 
perpendicular (cx sin 77) and longitudinal (cx cos 77) terms 
in (^1^) add up independently. 

In the environment-dominated regime, AE <^ ak^T, 
the qubit's Hamiltonian is too weak to fix the basis, while 
the coupling to the bath becomes the dominant part of 
the total Hamiltonian. Therefore one should discuss the 
problem in the eigenbasis of the observable Cz to which 
the bath is coupled. The spin can tunnel incoherently 
between the two eigenstates of az- To find the tunnel- 
i ng ra te one can again use the can onica l transformation 
( 4.13 ) leading to the Hamiltonian ( 1.14 ). In the Golden 
Rule approximation one obtains (Leggett et al. 1987) the 
following relaxation rate (for AE ^ akBT): 
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(4.18) 



where P{..) is the Fourier transform of P{t). To find the 
dephasing rate we use again Eq. ( [4.15 ) and obtain 
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(4.19) 



In this regime the dephasing is much faster than the re- 
laxation. Moreover, we observe that t~1^^ cx {akBT)~^ 
while oc akBT. This implies that the faster is the 
dephasing, the slower is the relaxation. Such a behavior 
is an indication of the Zeno (watchdog) effect (Harris & 
Stodolsky 1982): the environment frequently "observes" 
the state of the spin, thus preventing it from tunneling. 



C. Several fluctuating fields and many qubits 

Next we consider the general case of a qubit coupled 
to several fluctuation sources (baths) via different spin 
components and below we discuss a many-qubit system. 
It is described by the following generalization of the spin- 
boson model: 



(4.20) 



where the index i labels the different (Ohmic) baths 
and the unit vectors determine the spin components 
to which the baths are coupled. The Hamiltonian- 
dominated regime is realized when AE 3> X^i onkBT, 
where Ui correspond to the bath i. In this case one should 
divide the fluctuations of all baths into transverse and 



longitudinal ones, as in Eq. (4.8), with each bath being 
characterized by the angle 77^ between rii and the field di- 
rection (cos77i — Bni/\B\). The transverse fluctuations 
will add up to the relaxation rate as 
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while the longitudinal fluctuations lead to the dephasing 
rate 



' 9 relax ' / j 



na, cos^ 77, . (4.22) 



In the simplest environment-dominated situation one 
of the baths {i — io) is much stronger than all others, 
Q^io ^ '^i^io '^i^ ^^'^ satisfies ai„k^T ^ AE. Then, the 
dephasing is described in the eigenbasis of crrii^ corre- 
sponding to this bath. The rate is 



(4.23) 



The relaxation rate in this basis may be estimated using 
the Golden Rule: 



r-^ + > 27ra, sin^ Xi 



knT 



(4.24) 



where cosxi = (n^n j,,). This rate is smaller than the 
dephasing rate ( 4.23| ). Note that the relaxation rates 
due to the other than the strongest baths do not show 
the Zeno effect. 

In the most complicated case of the environment- 
dominated regime with several baths coupled to different 
spin components and characterized by constants ai of the 
same order it is difficult to make quantitative predictions. 
However, we expect that the time scale for a quantum 
state destruction (either dephasing or relaxation) will be 
longer than (^. liraik-QT /h)~ . 

So far we were concerned with dissipative effects in a 
single qubit. However, any but the simplest applications 
of quantum state engineering make use of many coupled 
qubits and entangled states, i.e., states whose properties 
can not be reduced to the single-bit ones. Therefore the 
question arises how dissipation affects multi-qubit sys- 
tems and entangled states. 

As a first step we analyze the effect of dissipation on an 
A/'-qubit system during an idle period when the single-bit 
terms and two-bit interactions are switched off, Tictii = 0. 
We assume that each qubit is coupled to an independent 
oscillator bath, with no correlations between the baths, 
and we choose a basis where this coupling is diagonal, 
i.e., the bath is coupled to the (Tx-component. In this 
case the time evolution operator, governing the density 
matrix, Pvi..MN:i-Li---tiN{^)^ factorizes. We perform for each 
qubit i the unitary rotation (4.13) with the result: 



N 



i=l 



(4.25) 



Averaging over the baths yields in the long-time limit the 
following time dependence of the density matrix: 

{t) (X n (-*/^^) ■ (4.26) 

Here the product is over all qubits i which have off- 
diagonal entries in the density matrix. This form shows 
that the dissipation has the strongest effect on those en- 
tries of p which are off-diagonal with respect to each 
qubit. For instance, the dephasing rate of po...O;i...i(0 
is = 1/t^- It scales linearly with the number of 
qubits. 



The result (4.26) applies independent of the initial 
state, i.e., equally for product states or entangled states 
of the multi-qubit system. However, it is valid only if the 
controlled parts of the Hamiltonian are switched off. The 
question how dissipation influences the dynamics of en- 
tangled states in general situations, e.g., during logic op- 
erations when the many-qubit Hamiltonian is non-zero, 
remains open. Further work is needed to analyze this 
interesting and important problem. 



D. Dephasing in charge qubits 

We now turn to the specific case of a Josephson charge 
qubit coupled to the environment and determine how the 
dephasing and relaxation rates depend on system param- 
eters. The system is sensitive to various electro-magnetic 
fluctuations in the external circuit and the substrate, as 
well as to background charge fluctuations. We first esti- 
mate the effect of fluctuations originating from the circuit 
of the voltage sources. In Fig. |l^ the equivalent circuit 
of a qubit coupled to an impedance Z(ijj) is shown. The 
latter has intrinsic voltage fluctuations with Johnson- 
Nyquist power spectrum. When embedded in the circuit 
shown in Fig. |l^ but with Ej = 0, the voltage fluctua- 
tions between the terminals of Z{u!) are characterized by 
the spectrum 

f huj \ 



(SVSV)^ = Re{Zt{uj)}hLo coth 



V 2fcBT / 



(4.27) 



Here Zt{uj) = [iwCqb + -^"^(w)] is the total 
impedance between the terminals of Z(uj) and Cqb is the 
capacitance ( ^.15 ) of the qubit in the circuit. 
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FIG. 13. A qubit in an electro-magnetic environment charac- 
terized by the impedance Z{ij). 

Following Caldeira & Leggett (1983) we model the dis- 
sipative element Z{ll}) by a bath of harmonic oscillators 
described by the Hamiltonian Hb as in Eq. (1.4) (Shnir- 
man et al. 1997). The voltage fluctuations between the 
terminals of Z{uj) are represented by 6V — Xa^a, 
while the spectral function (4.6) has to be chosen as 
J{uj) — ujKe{Zt{Lo)} in order to reproduce the fluctua- 
tion spectrum (4.27). 

To derive the Hamiltonian wc introduce an auxiliary 
variable, the charge q on the gate capacitor, and add a 
term which couples q to the bath (a natural choice since 
the current through the impedance is q) . For simplicity of 
derivation it is convenient also to add a small inductance 
in series with the impedance Z to provide a 'mass' for the 
q-mode. Allowing for a time-dependent external voltage 
Vg(t) and integrating out q, we find the Hamiltonian: 
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[2en-CgV-g(t)] 
2(Cj + Cg) 



Ej cos 9 + 



2ma 



E 



(4.28) 



The voltage (t) in the last term (coupling to the bath) 
is relevant only if it depends on time and is usually 
dropped in the literature. Its role is to provide the RC 
time delay. This means that when Vg is changed the 
qubit feels the change only after the RC time. In the 
systems described here this delay is much shorter than 
the other relevant time scales and, thus, need not to be 
discussed further. 

To be specific we will concentrate in the following on 
the fiuctuations due to an Ohmic resistor Z{ll}) = Ry 
in the bias voltage circuit. In the two-state approx- 
imation, using the relations n = (1/2)(1 -I- tr^) and 
cos 6 = (l/2)f72;, we thus arrive at the Hamiltonian of 
the spin -boson mo del ( fl.l| ) with and B^- given by 
Eq. (2^) and ( |2.5| ), respectively. The dimensionless pa- 



rameter , characterizing the effect of fluctuations cou- 
pling to <Tz of the charge qubit, is given by 
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(4.29) 



The circuit resistance is compared to the quantum resis- 
tance i?K = h/e'^ « 25.8 kQ. Since the parameter af^ di- 
rectly relates the dephasing rate to typical energy scales, 
its inverse determines the number of coherent single- 
qubit manipulations which can be performed within the 
dephasing time. From Eq. ( 4.29 ) we see that in order 
to keep the dissipative effects of external voltage fiuctua- 
tions weak one has to use a voltage source with low resis- 
tance and choose the gate capacitance Cg « Cqb <C Cj 



as low as possible. The latter screens out the voltage 
fiuctuations, at the expense that one has to apply larger 
gate voltages for the manipulations. At high frequencies 
discussed here the typical impedance of the voltage cir- 
cuit is Ry « 50 O-jand one obtains af « 10-^{Cg/Cj)^ . 
Taking the ratio Cg/Cj = 10^^ one can reach the dis- 
sipation as weak as af^ « 10^^, allowing in principle for 
10^ coherent single-bit manipulations. 

The fluctuations of the externally controlled flux 
through the SQUID loop of a charge qubit with tunable 
Josephson coupling (see Fig. ^ also destroy the phase 
coherence. As can be seen from Eq. ( p. 10 ) these fluctu- 
ations couple to cos Q (X ax- Analogously to the voltage 
fluctuations considered above, their strength can be ex- 
pressed by the effective impedance of the current circuit 
which supplies the flux. For an estimate we take this 
impedance Rj to be purely real. At typical high frequen- 
cies of the qubit 's operation it is of order of the vacuum 
impedance, Rj ^ 100 O. In terms of this resistance the 
effect of fluctuations in the flux is characterized by the 
parameter 
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(4.30) 



where the flu x-con trolled Josephson coupling i?j(<l'x) is 
given by Eq.( 2.11 ). The effect is weak for low values 
of the mutual inductance. For M k, 0.01 — 0.1 nH and 
£;5 « 0.1 K we obtain af « 10"^ - 10"*^. The dephas- 
ing and relaxation times, if only flux fluctuations need 
t o be considered, are thus given by Eqs. ( 4.11 ), ( 1.12 ), 
( 4.30 ), but with a substitution tan 77 = Bz/B.^ since the 
noise terms couple to cTx- When both the gate voltage 
and the flux fluctuate a multi-bath situation described 



by Eq. ( [4.20D is reahzed. 

For typical parameters, e.g., for those of the experi- 
ments of Nakamura et al. (1999), one can estimate the 
dephasing time associated with the noise in the external 
circuit to be of the order of 100 ns. These experiments 
allow a direct probing of the phase coherence. Coher- 
ent oscillations have been observed for 5 ns. Hence Xhe 
theoretical estimate appears to be of the right order U. 



^"Nakamura et al. (1999) reached an even smaller ratio 
for the qubit, but the probe circuit introduced a high stray 
capacitance. 

^^In the experiments of Nakamura et al. (1999) actually 
much of the dephasing can be attributed to the measurement 
device, a dissipative tunnel junction which was coupled per- 
manently to the qubit. Its tunneling resistance was optimized 
to be large enough not to destroy the qubits quantum coher- 
ence completely, but low enough to allow for a measurable 
current. Single-electron tunneling processes, occurring on a 
time scale of the order of 10 ns, destroy the state of the qubit 
(escape out of the two-state Hilbert space) thus putting an 
upper limit on the time when coherent time evolution can be 
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Another important source of decoherence in charge 
qubits are fluctuations of the background charge. It 
was found experimentally that they lead to 1/f noise 
at low frequencies. Typically, their contribution to fluc- 
tuations of the effective gate voltage, Qg = CgVg, is of 
order 5'Q(tj) — lO^^e^/w. These fluctuations Hmit the 
time of coherent evolution and can hardly be improved 
by present-day experimental techniques 113. Fortunately, 
for the ideas of quantum state engineering, the configu- 
rations of background charges change rarely, on a typi- 
cal slow relaxation time scale which can reach minutes 
or even hours. Thus during each cycle of manipula- 
tions these charges provide random static gate voltages, 
which do not destroy, e.g., coherent oscillations. Their 
effect can be suppressed using refocusing techniques. On 
the other hand, when averaging over many experimen- 
tal runs (as done by Nakamura et al. (1999)), then an 
ensemble-average is formed and the dephasing rate can 
be extracted. 

When the qubit is coupled to a measurement device, 
which necessarily implies a coupling to a macroscopic 
variable with dissipative dynamics, the feedback intro- 
duces fluctuations and causes dephasing. We will discuss 
this explicitly in the next section for the case where a 
Josephson charge qubit is coupled to a dissipative single- 
electron transistor. We find that the shot noise of the 
tunneling current in the SET introduces dephasing and 
relaxation processes. 



E. Dephasing in flux qubits 

Flux qubits have the advantage that they are practi- 
cally insensitive to the background charge fluctuations 
and are advertised for this reason (Mooij et al. 1999). 
Still their phase coherence can be destroyed by a num- 
ber of effects as well. Some sources of dissipation for 
flux qubits have been discussed by Tian et al. (2000) 
and estimates have been provided for the parameters of 
the circuits of Mooij et al. (1999) and Orlando et al. 
(1999). This includes the effect of background charge 
fluctuations (r^ « 0.1 s) as well as quasiparticle tun- 
neling in the superconductor with a non-vanishing sub- 
gap conductance (r^ ~ 1 ms). The effect of nuclear 
spins in the substrate producing fluctuating magnetic 
fields is similar to the effect of background charges on 
charge qubits. While the static random magnetic field 
may induce substantial changes of the qubit frequencies 



observed. For a more detailed discussion of the experiment 
and the measurement process we refer to the article by Choi 
et al. (2000). 

'^■^See, however, recent work of Krupenin et al. (2000) where 
the l//-noise was suppressed by fabricating a metallic island 
on top of an electrode instead of placing it on the substrate. 



of order i5f2nuci ~ 30 MHz (which can be suppressed, in 
principle, by using refocusing pulses), they cause no de- 
phasing for the qubits until a typical nuclear spin relax- 
ation time Ti, which can reach minutes. Other sources 
of dephasing studied by Tian et al. (2000) include the 
electro-magnetic radiation (r^p « 10'' s), much weaker 
than in typical rf-SQUlD designs, and unwanted dipole- 
dipole magnetic couplings between qubits, which for the 
inter-qubit distance of 10 /xm produces substantial effects 
after a relatively short time r « 0.2 ms. Variations of the 
design were suggested to reduce the latter effect. 

An important source of dissipation of flux qubits are 
the fluctuations in the external circuit which supplies 
fluxes through the loops. They can be analyzed along the 
same lines as presented above. Since they couple to CTz, 
the relevant parameter is a^. It is fixed by the impedance 
Zj of the current source in the input loop providing the 
flux bias, and the mutual inductance of the input and the 
qubit's loop, M: 



4 ^ 



(4.31) 



The numerical prcfactor in Eq. ( 4.31 ) is a = 6(/3l — l)/7r^ 
for the rf-SQUID and a — b'^y/Ab'^ — for the design 
of Mooij et al. (1999) with 6 = E.s/Ej being the ratio 
of critical currents for the junctions in the loop (see Sec- 
tion Ill.A). The dephasing is slow for small loops and 
junctions with low critical currents. Indeed, the argu- 
ment in the bracket in (4.31) can be represented as the 
product of the 'screening ratio' Al/ L and the quantity 
= An'^LEj/^l. While the latter should be larger 
than one for the rf-SQUID, it can be much smaller for 
the design of Mooij et al. (1999), resulting in a slower 
dephasing. 

The impedance Zj(lu) and the bath spectrum J{uj) = 
hujTie ZY^{lu) are frequency-dependent. They can lead 
to resonances, e.g., if the current source is attached to 
the qubit via low-loss lines. Care has to be taken in ex- 
periments to avoid these resonances. In addition, the 
dephasing and relaxation times have to be esti mate d for 
this situation. The analysis of the Hamiltonian ( [4.8| ) with 
general bath spectrum shows that they can be defined 
as the times when the quantities cos^ 77 ((J^ X(t)c?t)^) 
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reach values of or- 



and sin ri{\J^X{T)e^ 
der one, respectively. The first quantity can be ex- 
pressed as t J duj{X^)6{uj) involving the function 6{cj) = 
2sm^{iL>t/2)/{Trtuj^) which is peaked at w = and has 
the width Recalling that {X^) is related to the bath 
spectrum ([4.5| ), one finds that only the low- frequency 
part of the impedance Zi{uj), with uj < t^^, determines 
the dephasing rate. On the other hand, the relaxation 
rate depends on the values of Zj{uj) at frequencies in the 
vicinity of uj — AE/h in a range of width t~^^^. 

The choice of a high dc resistance of the remote cur- 
rent source strongly suppresses the fluctuations at low 
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frequencies. Using Eqs. (4.31, 4.12) we can estimate 
the dephasing rate for the parameters of Mooij et al. 
(1999). Assuming w 0) - 1 Mfi and M « L we 

find ~ 10~^, which implies a negligible dephasing. 
At higher frequencies of order Ai?, even if resonances 
are suppressed, Zi{tJ ~ AE) is of order 100 il, leading 
to ~ 10~^. This determines the relaxation rate via 



Eq. (4.11). For Mooij et al.'s (1999) parameters we es- 
timate a relaxation time of ~ 5 /iS at the degeneracy 
point. 

The effect of fluctuations of the crj;-term in the Hamil- 
tonian (in the flux circuit of the dc-SQUID-loop which 
controls the Josephson coupling) can be described in a 
similar way (Makhlin et al. 1999). The effect of these 
fluctuations is relatively weak for the operation regimes 
discussed by Mooij et al. (1999). 

While a further analysis of the dephasing effects in flux 
qubits may be needed, the above-mentioned estimates of 
suggest that the observation of coherent flux oscilla- 
tions is feasible in the near future. However, an open 
problem remains the fact that the observation requires 
an efficient quantum detector, e.g., a quantum magneto- 
meter. We discuss this issue in the next section. 



V. THE QUANTUM MEASUREMENT PROCESS 
A. General concept of quantum measurements 

Quantum state engineering requires controlled quan- 
tum manipulations but also quantum measurement pro- 
cesses. They are needed, e.g., at the end of a computa- 
tion to read out the final results, or even in the course 
of the computation for the purpose of error correction. 
The problem of quantum measurement has always at- 
tracted considerable attention, and it still stirs contro- 
versy. In most of the literature on quantum information 
theory the measurement process is expressed simply as 
a "wave function collapse", i.e., as a non-unitary projec- 
tion, which reduces the quantum state of a qubit to one 
of the possible eigenstates of the observed quantity with 
state-dependent probabilities. 

On the other hand, in reality any measurement is per- 
formed by a device which itself is a physical system, suit- 
ably coupled to the measured quantum system and with 
a macroscopic read-out variable. This is accounted for in 
the approaches described below, which are based on the 
concepts of the dissipative quantum mechanics (Leggett 
et al. 1987, Weiss 1999, Zurek 1991). The qubit and the 
measuring device are described as coupled quantum sys- 
tems. Initially there exist no correlations between the 
qubit and the meter but, due to the coupling, such corre- 
lations (entanglement) emerge in time. This is precisely 
the concept of measurement of von Neumann (1955). 
Furthermore, the meter is actually a dissipative system, 
coupled to the environment. The dissipation eventually 



reduces the entanglement between the qubit and the me- 
ter to classical correlations between them. This is exactly 
what is needed for the measurement. 

Several groups have studied problems related to quan- 
tum measurement processes in mesoscopic systems. 
Aleiner et al. (1997), Levinson (1997), and Gurvitz 
(1997) investigated the effect of a dissipative conductor, 
whose conductance depends on the state of quantum sys- 
tem, onto this system itself. Such a configuration is re- 
alized in recent experiment (Buks et al. 1998, Sprinzak 
et al. 1999), where a quantum dot is embedded in one of 
the paths of a 'which path' interferometer. The flow of a 
dissipative current through a nearby quantum point con- 
tact, with conductance which depends on the charge on 
the dot, amounts to a 'measurement' of the path chosen 
by the electrons in the interferometer. And indeed the 
flow of a current through the point contact leads to an ob- 
servable reduction of the flux-dependence of the current 
through the interferometer. On the other hand, no real 
measurement was performed. The experiments merely 
demonstrated that the amount of dephasing can be con- 
trolled by a dissipative current through the QPC, but 
they did not provide information about the path chosen 
by each individual electron. 

In what follows we will discuss systems in which one 
is not only able to study the dephasing but also, e.g., by 
measuring a dissipative current, to extract information 
about the quantum state of the qubit. As an explicit ex- 
ample we investigate a single-electron transistor (SET) in 
the sequential tunneling regime coupled capacitively to a 
charge qubit (Shnirman & Schon 1998, Schoelkopf et al. 
1998). Our purpose is not philosophical but rather prac- 
tical. We do not search for an ideal measurement device, 
instead we describe the properties of a realistic system, 
known to work in the classical regime, and investigate 
whether it can serve as quantum measurement device. At 
the same time we keep in mind that we should come as 
close as possible to the projective measurement picture 
assumed in quantum algorithms. With the same goal, 
namely, to establish their potential use as quantum detec- 
tors, various other mesoscopic devices have been studied 
recently. These include a quantum point contact (QPC) 
coupled to a quantum dot (Gurvitz 1997, Korotkov 1999), 
a single-electron transistor in the co-tunneling regime 
(Averin 2000 a, Maassen van den Brink 2000), or a su- 
perconducting SET (SSET) (Averin 2000 &, Cottet et al. 
2000). For a flux qubit a dc-SQUlD, or suitable modifi- 
cations of it, can serve as a quantum detector. 

Technically the measurement process is described by 
the time evolution of the reduced density matrix of the 
coupled system of qubit and meter. To analyze it we first 
derive a Bloch-master equation for the time-evolution of 
the system of a Josephson charge qubit and a dissipa- 
tive single-electron transistor (SET) coupled to it. This 
derivation demonstrates how the unitary time evolution 
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of system plus detector can lead to a quantum measure- 
ment. It also allows us to follow the dynamics of the 
system and detector and to analyze the mutual influence 
of their variables. During the measurement the qubit 
looses its phase coherence on a short dephasing time scale 
This means that the off-diagonal elements of the 
qubit's density matrix (in a preferred basis, which de- 
pends on relative strengths of the coupling constants) 
vanish, while the diagonal ones still remain unchanged. 
At the same time the information about the initial state 
of the qubit is transferred to the macroscopic state of 
the detector (the current in the SET). Under suitable 
conditions, after another time scale Tmeas, this informa- 
tion can be read out. Finally, on a longer, mixing (or 
relaxation) time scale Tmix the detector acts back onto 
the qubit and destroys the information about the initial 
state. The diagonal entries of the density matrix tend to 
their stationary values. (These are either determined by 
the detector or they are thermally distributed, depending 
on the relative strength of the measurement device and 
residual interactions with the bath.) One has to chose 
parameters such that this back-action does not change 
the occupation probabilities of the qubit's state before 
the information is actually read out, T,„ix ^ T,„eas > t^. 

The different times scales characterizing the measure- 
ment process by a SET, show up also in equilibrium prop- 
erties, e.g., in the noise spectrum of the fluctuations in 
the SET. This can be extracted from the time evolution 
of the coupled density matrix as well. 

Ideally, the meter is coupled to the qubit only during 
the measurement. In practice, however, this option is 
hard to realize for mesoscopic devices. Rather the meter 
and the qubit are coupled permanently but the former 
is kept in a non-dissipative state. To perform a mea- 
surement, the meter is switched to the dissipative regime 
which, as we mentioned above, is an important require- 
ment for the effective quantum measurement. For in- 
stance, in a SET, due to the Coulomb blockade phenom- 
ena, no dissipative currents are flowing in the meter as 
long as the transport voltage is switched off. Applying 
a voltage bias above the Coulomb threshold induces a 
dissipative current through the SET, which leads to de- 
phasing and, at the same time, provides the macroscopic 
read-out variable. 

Various ways of operation can be used, depending on 
details of a particular setting. For instance, one can 
switch the meter abruptly into the dissipative regime and 
monitor the response of dissipative currents to the qubit. 



Another possibility is to change the bias gradually un- 
til the system switches into the dissipative regime. The 
value of the bias at which the switching occurs provides 
information about the qubit's state. Conceptually, these 
techniques are similar. For definiteness, below we discuss 
the former operation strategy. 

B. Single-electron transistor as a quantum electrometer 

Since the relevant quantum degree of freedom of a 
Josephson charge qubit is the charge of its island, the 
natural choice of measurement device is a single-electron 
transistor (SET). This system is shown in Fig. ^ The 
left hand part is the qubit, with state characterized by 
the number of extra Cooper pairs n on the island, and 
controlled by its gate voltage V^''''. The right hand part 
shows a normal island between two normal leads, which 
form the SET. Its charging state is characterized by the 
number of extra single-electron charges N on the middle 
island. It is controlled by gate and transport voltages, 
V^'^^ and Vtr, and further, due to the capacitive cou- 
pling to the qubit, by the state of the latter. A similar 
setup has been studied in the experiments of Bouchiat 
(1997) and Bouchiat et al. (1998), where it was used to 
demonstrate that the ground state of a single Cooper pair 
box is a coherent superposition of different charge states. 

During the quantum manipulations of the qubit the 
transport voltage Vtr across the SET transistor is kept 
zero and the gate voltage of the SET, Vg^^""", is chosen to 
tune the island away from degeneracy points. Therefore 
at low temperatures Coulomb blockade effects suppress 
exponentially a dissipative current flow in the system, 
and the transistor merely modifies the capacitances of 
the system lU. To perform a measurement one tunes the 
SET by V^^'^ to the vicinity of its degeneracy point and 
applies a small transport voltage Vtr- The resulting nor- 
mal current through the transistor depends on the charge 
configuration of the qubit, since different charge states in- 
duce different voltages on the island of the SET. While 
these properties are well established as operation prin- 
ciple of a SET as electro-meter in the classical limit, it 
remains to be demonstrated that they also allow resolv- 
ing different quantum states of the qubit. For this pur- 
pose we have to discuss various noise factors, including 
the shot noise associated with the tunneling current and 
the measurement induced transitions between the states 
of the qubit. They can be accounted for by analyzing 
the time evolution of the density matrix of the combined 
system. 



^ In this section denotes the dephasing time during a 
measurement. It is usually much shorter than the dephas- 
ing time during the controlled manipulations discussed in the 
previous sections, from the context it should be clear which 
situation we refer to. 



^ More precisely, the leading contribution are co-tunneling 
processes, which are weak in high-resistance junctions. 
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FIG. 14. The circuit consisting of a qubit plus a single- 
electron transistor used as a measuring device. 

The Hamiltonian of the combined system consists of 
the parts describing the qubit, the SET and the interac- 
tion between them: 



W — Wctil + WsET + Wii 



(5.1) 



Except for a redefinition and renormalization of parame- 
ters, the qubit's part is identical to that of the bare qubit 



with = 4£'c,qb(l - 2ng,qb) and = Ej{<^^). For a 
detailed derivation and precise definition of the parame- 
ters see Appendix |^. Here it is sufficient to know that 
the qubit's Hamiltonian can be controlled by gate volt- 
ages an d th e flux through the SQUID. After diagonaliza- 
tion of (|2!|), Hctri = -(l/2)A£;(r;) p^, the characteristic 



energy scale of the qubit, AE{r]) — \/ B'^ 
apparent. 

The Hamiltonian of the SET reads 

'^SET — -Ec,SEt(A^ — A^g.SEx)^ 

+ Hh + + Hi+ TLt , 



B^, becomes 



(5.2) 



where the transistor's charging energy is given by (C2) 
and the gate charge iVg^sET = — eV^/2i?c,SET, defined 
in Eq. (p4D, can be controlled by Vg^^""". The three terms 
T^L, Ti-R, and Tii describe microscopic degrees of freedom 
of noninteracting electrons in the two leads and the mid- 
dle island of the SET transistor: 



(r = L,R,I) . 



(5.3) 



The index r = L, R, I labels the electrodes "left" , "right" 
(viewed from a suitable angle) and island, a labels trans- 
verse channels including the spin, while k refers to the 



wave vector within one channel. Note that similar terms 
should have been written for the electrode and island 
of the qubit; however, for this superconducting non- 
dissipative element the microscopic degrees of freedom 
can be integrated out (Ambegaokar et al. 1982, Schon 
& Zaikin 1990), resulting in the "macroscopic" quan- 
tum description presented in Sections |l| and HI. In 
this limit the tunneling terms reduce to the Josephson 
coupling Hj = —Ej cos 9, expressed in a collective vari- 
able describing the coherent transfer of Cooper pairs. 



n) 



1 



The normal-electron tunneling in the SET transistor is 
described by the standard tunneling Hamiltonian, which 
couples the microscopic degrees of freedom: 



kk' (7 

+ ^ k' k" a'^k" a^k' 



iff) 



(5.4) 



k' k" (7 



To make the charge transfer explicit, (5.4) displays two 
"macroscopic" operators, e^"^ and e^^\ The first one 
describes changes of the charge on the transistor island 
due to the tunneling: e"^\N) = |A^+1). If the total num- 
ber of electrons on the island is large it may be treated 
as an independent degree of freedom. We further include 
the operator e^"^"^ ^ which acts on m, the number of elec- 
trons which have tunneled through the SET transistor, 
e^^\m) ~ 1 771 + 1). Since the chemical potential of the 
right lead is controlled, 777 docs not appear in any charg- 
ing part of the Hamiltonian. However, we have to keep 
track of it, since it is the measured quantity, related to 
the current through the SET transistor. 

Finally, Hint — Ei^tN{2n — 1) describes the capacitive 
interaction between the charge on the qubit's island and 
the SET island. In detail it originates from the mixed 
term in Eq. (CI), where E\at is given by Eq. (|C2|). In the 
2-state approximation for the qubits, 77 = (1 + Gz)/'^, we 
obtain 



H 



int 



N5n 



int 



NE^ 



int CT z 



(5.5) 



The operator <57iint, introduced here for later conve- 
nience, is the part of the interaction Hamiltonian which 
acts in the qubit's Hilbert space. 



C. Density matrix and description of measurement 

The total system composed of qubit and SET is de- 
scribed by a total density matrix g(t). We can reduce it, 
by tracing out the microscopic electron states of the left 
and right leads and of the island, to 



m - TrL.Rj{e(i)} . 



(5.6) 



This reduced density matrix p{i, i'\ N, N'; 777, 777') is still a 
matrix in the indices i, which label the quantum states of 
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the qubit |0) or |1), in N, and in m. In the following we 
will assume that initially, as a result of previous quantum 
manipulations, the qubit is prepared in some quantum 
state and it is disentangled from the SET, i.e., the initial 
density matrix of the whole system may be written as a 
product po = (g) Pq^^- At time t — we switch on 
a transport voltage of the SET and follow the resulting 
time evolution of the density matrix of the whole sys- 
tem. For specific questions we may further reduce the 
total density matrix in two ways, either one providing 
complementary information about the measurement pro- 
cess. 

The first, common procedure (Gurvitz 1997) is to trace 
over N and m. This yields a reduced density matrix of 
the qubit pij = ^p{i,i]N,N;ra,m). At t = one 

has Pij = {p^)ij- Depending on the relation among the 
energy scales and coupling strengths (see Section IV) a 
preferred basis may exist in which the dynamics of the di- 
agonal and the off-diagonal elements of pij decouple. We 
study how fast the off-diagonal elements (in that special 
basis) vanish after the SET is switched to the dissipative 
state, i.e., we determine the rate of dephasing induced by 
the SET. And we determine how fast the diagonal ele- 
ments change their values. This process we call 'mixing'. 

The rates of 'dephasing' and 'mixing' refer to the quan- 
tum properties of the measured system in the presence 
of the measurement device. These quantities have also 
been analyzed in Refs. (Levinson 1997, Aleiner et al. 
1997, Gurvitz 1997, Buks et al. 1998). They do not tell 
us, however, anything about the quantity measured in 
an experiment, namely, the current flowing through the 
SET. Therefore the second way to reduce the density ma- 
trix is important as well. By tracing the density matrix 
over the qubit's variables and the state N of the island. 



P(m, t) = ^ p{t, I- N, N; m, m){t) , 

i,N 



(5.7) 



we obtain the probability distribution for the experimen- 
tally accessible number of electrons m which have tun- 
neled through the SET during time t. Its detailed analy- 
sis will be presented below. In order to provide a feeling 
we first present and describe some representative results: 
At t = no electrons have tunneled, so P{m, 0) = S,n.o- 
Then, as illustrated in Fig. ^ the peak of the distribu- 
tion moves to nonzero values of m and, simultaneously, 
widens due to shot noise. If the two states of the pre- 
ferred basis correspond to different tunneling currents, 
and hence different m-shifts, and if the mixing in this 
basis is sufficiently slow, then after some time the peak 
splits into two with weights corresponding to the initial 
values of the diagonal elements of p^^ in the preferred ba- 
sis (we will denote these values |ap and |6p = 1 — |ap). 
Provided that after sufficient separation of the two peaks 
their weights are still close to the original values, a good 
quantum measurement can be performed by measuring 



m. After a longer time, due to transitions between the 
states of the preferred basis (mixing) the two peaks trans- 
form into a broad plateau. Therefore there is an optimum 
time for the measurement, such that, on one hand, the 
two peaks are separate and, on the other hand, the mix- 
ing has not yet influenced the process. 




t[ns] 



FIG. 15. P{m,t), the probability that m electrons have tun- 
neled during time t (measured in nanoseconds). The initial 
amplitudes of the qubit's states are \a\ ^ = 0.75, |fe|^ = 0.25. 
E} = 0.1 K, the remaining parameters are given later in the 
text. 

Having introduced the relevant time scales, we can 
mention that the which-path interferometry (Buks et al. 
1998) can be thought of as a very short measurement 
process. Each particular electron ~ the observed quan- 
tum system ~ spends only a short time within the dot. 
This time of interaction with the meter may be shorter 
than the dephasing time. Therefore the coupling leads 
only to a slight suppression of the interference pattern. 



D. Master equation 

The time evolution of the density matrix leads to 
Bloch-type master equations with coherent terms. Exam- 
ples of this type have recently been analyzed in various 
contexts (Nazarov 1993, Schoeller & Schon 1994, Stoof 
& Nazarov 1996, Gurvitz & Prager 1996, Gurvitz 1997). 
Schoeller & Schon (1994) developed a diagrammatic tech- 
nique which provides a formally exact master equation 
for a SET as an expansion in the tunneling term Tix, 
while all other terms, including the charging energy, con- 
stitute the zeroth order Hamiltonian H.^. The time evo- 
lution of the reduced density matrix, p{t) = p(0)n(0,t) 
is expressed by a propagator n(t',t), which is expanded 
and displayed diagrammatically and finally summed in 
a way reminiscent of a Dyson equation. Examples are 
shown in Fig. |l^ and Appendix In contrast to equilib- 
rium many-body expansions, since the time dependence 
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of the density matrix is described by a forward and a 
backward time-evolution operator, there are two prop- 
agators, which are represented by two horizontal lines 
(Keldysh contour). The two bare lines describe the co- 
herent time evolution of the system. They are coupled 
due to the tunneling in the SET. The sum of all dis- 
tinct transitions defines a 'self-energy' diagram S. Below 
we will present the rules for calculating S and present a 
suitable approximate form. The Dyson equation is equiv- 
alent to a Bloch-master equation for the density matrix, 
which reads 



dt'p{t')^{t - t') 



(5.8) 



n 



n 



FIG. 16. The Dyson-type equation governing the time evolu- 
tion of the density matrix. It is equivalent to the generalized 
master equation (5.8). The 'self-energy' diagrams E describe 



the transitions due to tunneling in the SET transistor. 

In principle, the density matrix p{i, i' \ N, N'; m, m') = 

Pi 'n 'm ^ matrix in all three indices j, N, and m, 
and the (generalized) transition rates due to single- 
electron tunneling processes (in general of arbitrary or- 
der), '"^ ^V^, '™ (t—f), connect these diagonal and 
off-diagonal states. But a closed set of equations for the 
time evolution of the system can be derived (Schoeller & 
Schon 1994) which involves only the elements of p diag- 
onal in N. The same is true for the matrix structure in 
771. I.e., we need to consider only the following elements 
of the density matrix p^{m) = p*'^'™. Accordingly, of 
all the transition rates we need to calculate only the cor- 

,N' ,m.' ^i,N 



responding matrix elements S -, 



(At). 



If the temperature is low and the applied transport 
voltage not too high, the leading tunneling processes in 
the SET are sequential transitions between two adjacent 
charge states, say, = and A*" = 1. We concentrate 
here on this case (to avoid confusion with the states of 
the qubit we continue using the notation A^ and A^ -I- 1). 
The transition rates can be calculated diagrammatically 
in the framework of the real-time Keldysh contour tech- 
nique. The derivation is presented in Appendix A 
simplification arises, since the coherent part of the time 
evolution of the density matrix, which evolves on the 
time scale set by the qubits energies, is slow compared 
to the 'tunneling time', given by the inverse of the en- 
ergy transfer in the tunneling process (see arguments 
in square brackets in Eqs. 5.11, ^.12| below). As a re- 
sult the self-energy effectively reduces to a delta-function, 
T,{At) oc d{At), and the Bloch-master equation (|5.8|) re- 



duces to a Markovian dynamics. 

The resulting master equation is translationally invari- 
ant in m-space. Hence a Fourier transformation is appro- 
priate, /5^(fc) = ^ .-^km;,^3 



pfj{m). As a result Eq. ( p8| ) 



factorizes in fc-space and we get a finite rank (8 x 8) sys- 
tem of equations for each value of k. This system may be 
presented in a compact form if we combine the eight com- 
ponents of the density matrix into a pair {pN, PN+i) of 
the 2 X 2-matrices p^{k), corresponding to A^ and A^ + 1: 



4 

at V pn^ 



Pn 



i[Wctrl,/57v] \ 
i[Wctrl + SHiut, Pn+i] J 



El 



(5.9) 



The operator STi, 



int 



EintO-z — E-^ntP'^ 



E^^Px was in- 



troduced in Eq. (5.?;). The tunneling rates in the left and 
right junctions are represented by operators Fl and Fr, 
acting on the qubit 's density matrix: 



ErPW+i = r^PN+l — TTQIR [STiint^PN^ 



(5.10) 



Here q;l/r = i?K/(47r^i?j!;^j^) is the tunneling conduc- 
tance of the left /right junction, measured in units of the 
resistance quantum i?K — h/e^. The tunneling rates in 
the junctions are determined by the potentials Pl/r of 
the leads and the induced charge A^g,SET on the SET's 
island: 

Fl = 2^aL[pL - (1 - 2A^g^sET)Sc,SET] (5.11) 
Fr = 2^aR[(l - 2Arg,sET)Sc,SET - Mr] (5-12) 
F = FlFr/(Fl+Fr). (5.13) 

They combine into the parameter F which gives the con- 
ductance of the SET in the classical (e.g., high voltage) 
regime. As a result of the last terms in the Eqs. ( 5.10| ) 
the effective rates and thus the current in the SET are 
sensitive to the state of the qubit, which makes the mea- 
surement possible. 
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FIG. 17. P{m,t), the probability that m electrons have tun- 
neled during time t. The parameters are the same as in Fig. 
except that Ej = 0.25 K. 

To provide a feeling for the types of solutions we 
present here numerical solutions of Eq. (5.9) for the fol- 
lowing system parameters: Bz = 2 K, Emt = 0.25 K, 
aL = aR = 0.03, Tl = 1.8 K, and Tr, = 7.8 K. We plot 
the resuhs for £'j = 0.1 K (Fig. ^ and for Ej = 0.25 K 
(Fig. |l^. We see that for the smaller value of £'j (the 
eigenbasis of Tictri is closer to the charge basis) the proba- 
bility distribution P(m, t) develops a two-peak structure. 
The weights of the peaks are equal to the initial values of 
the diagonal elements of the qubit's density matrix. For 
the larger value of Ej the two-peak structure can also 
be seen but the valley between the peaks is filled. This 
indicates that the mixing transitions take place on the 
same time scale as that of peak separation, and no good 
measurement can be performed. In Fig. |l8| the probabil- 
ity distribution P(m,t) is plotted for longer times. (In 
order to cover many decades of the time it is necessary to 
rescale the m axis as well.) The parameters are the same 
as in Fig. ^ {Ej = 0.1 K). The figure displays clearly 
the measurement stage and the third stage on the mix- 
ing time scale. The valley between the two peaks fills 
up and a single, broad (note the rescaled m axis) peak 
develops. 



As one can see from (5.9) the Hamiltonian of the qubit 
switches between Tictri and Tictri + SHint after each tun- 
neling event in the SET. This leads, in general, to a com- 
plicated dynamics. Above we showed repres ent ative nu- 
merical results. Next we will analyze Eq. (x£) pertur- 
batively which provides insight into the physics of the 
measurement process. 



P(m,t) 




me as 



FIG. 18. P{m, t) obtained in the two-mode approximation 



( 5.33), plotted versus t (on a logarithmic scale) and m/t. The 
initial probabilities of the qubit's states are \a\^ = 0.75, \b\^ = 
0.25. 



In general, the two Hamiltonians Tictri and 7ictri+<^7iint 
are not close and their respective eigenbases may be quite 
different. In this situation a better choice for the qubit's 
Hamiltonian would be the "average" Hamiltonian 



Ha 



n 



Ctrl 



{N)6H 



int 



(5.14) 



where {N) = rL/(rL + Fr) is the steady state av- 
erage of N . (The left-right asymmetry follows from 
the specific choice of the sequential tunneling regime: 
TV = iV = 1 in the fe/t junction and iV = 1 ^ TV = 
in the right junction.) Relevant energy scales are, then, 
the level splitting of the "average" Hamiltonian, 



Ha 



(5.15) 



the capacitive coupling energy Ei^t, and the bare tunnel- 
ing rates F^/^. For definiteness we consider the case 



Eint « Tl 



F, 



(5.16) 



I.e., we exclude the regimes of too strong coupling or too 
weak tunneling in the SET, which could be treated simi- 
larly to what is presented below. The main simplification 
of the regime (5.16) is that the spectrum of fiuctuations 
of the electron number N on the SET island is white 
in a wide enough range of frequencies, characterized by 
the zero frequency noise power S^ii^J ^ 0) = 2{SN^^q). 
When the SET is detached from the qubit and switched 
to the dissipative regime, this spectrum is given by 
SnIco ^ 0) = 4r/(FL + Fr)2. The related back-action 
noise randomizes the relative phase between the charge 
states of the qubit. In the absence of other sources for 
dynamics the rate of this process is 



'ipO 



(Fl + TrF 



(5.17) 



It is an important sca le for the distinction between dif- 
ferent sub-regimes of (l5J6|) , which will be described in 
the following. 



E. Hamiltonian-dominated regime 

Let us first consider the case where the qubit's "av- 
erage" Hamiltonian dominates over the dephasing due 



to the back-action of the SET, AE^v > t^q ■ Then, 
a perturbative analysis in the eigenbasis of Tiav is ap- 
propriate, where the operator SH-mt can be rewritten 
as SHint = £^Lpz + E^tPx, where Ef^^ = i^int cos r^av , 



E^. 



Eint sin 77av , and sin ry, 



where E"^^ = 

„ = Sj/ASa 



In the fol- 



lowing we treat the off-diagonal part of STCint, i-c, Er^^p^, 



(which leads to mixing, see below) pcrturbatively. In 
zeroth order, i.e., for E^^ — 0, the time evolution of 
[p'^ik), p^j_-^{k)) with different are decoupled from 
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each other. For diagonal elements, i 
obtain 



or 1, we 



.00/11 

Pn 

.00/11 

Pn+1 



,0/1 



r' 



,0/1 



nO/l 



.00/11 

Pn 

.00/11 

Pn+1 



(5.18) 



where F 



0/1 



Fl ± and F^ 



0/1 



We obtain four eigenmodes (oc e'*'*), two for each el- 
ement (00 or 11). Most interesting are small values 
fc <C 1. In this case, two modes with eigenvalue 
A!!(fc) w — iF*fc — i/T*A:^ and eigenvectors V^{k) given 
by PN{k)/PN+i{k) ~ rk/rL' describe waves in m-space 
propagating with group velocity 



F ' EE 



(5.19) 



The wave-packets widen with time as ^j2pVH due to 
shot-noise effects. The so called Fano factors. 



r^(FL +Fi,,^)/(FL + Fy2, 



(5.20) 



determine the suppression of the shot noise in the sequen- 
tial tunneling regime. The second pair of eigenmodes 
decay fast, A"(fc) w — (r'L + '^\d- -'^^^ their eigenvec- 
tors we obtain /5]^(fc)//57v-i-i(^) ~ This fast decay 
means that after a few tunneling events, on the time scale 
(ri + rk)"\ detailed balance p%(}i) lp% ^^(k) « V^jVl 
is established. For larger values of /c (fc ~ 1) one can 
check that both p^}^{k) and p]^_|_i(fc) decay fast. 
For the off-diagonal elements we obtain 



, / .01/10 

± { Pn 

\ .01/10 

^* V Pn+1 



-ifc 



(5.21) 



± 



il^EN^ 



.01/10 

Pn 

.01/10 
^'jV+1 



where A£'a 



and Ai?/\ 



N = AE'av -t- 2{N)El[^^_, ana ^£/jv+i 
2(1 — {N))Ef^^_. Again, there are two pairs of eigenmodes 
with eigenvalues such that A^°(A;) = [A^}.^(— fc)]*. One pair 

-(Fl+Fr) decay 



of eigenmodes with ReX'^^^^{k ^ 1) 



fast. For the second pair 
where 



ATE!' : II 2 

(1 L + 1 Rj^ 



T^a COS^ 77av 



(5.22) 



Thus after detailed balance is established, the AL" and A°^ 
modes describe coherent oscillations of the off-diagonal 
matrix elements with frequency of order Ai^av and decay 
rate r^^. For larger values of k the decay times are of 
the same order or shorter than those for fc « 0. 

The fast decaying diagonal modes A*|°^^^ do not con- 
tribute to P{m,t), since for these modes /5]^(A:) -f 



Plv+ii^) ~ 0. Thus there are only two modes contribut- 
ing, with eigenvalues A™(fc) and X^{k). Starting from 
the initial occupation probabilities |ap and |5p we obtain 



(5.23) 



where P{m,t) = J dk / {2tt) P {k , t) e'''"\ This form de- 
scribes the evolution of the distribution P(m, t) from the 
initial 5{m) at i = into two peaks. The peaks shift 
to positive m-values linear in t with velocities F° and 
F^, and they grow in widths as ^/2pTH. This time de- 
pendence implies that only after a certain time, which 
we denote as "measurement time" r^cas, the two peaks 
emerge from the broadened distribution. The associated 
rate is 



po_ri 



y^2f^+y/2rT^ 



(5.24) 



In the linear response regime, when F" and F^ are close, 
we obtain 



{All 
45/ 



(5.25) 



where Si is the zero frequency power of the shot noise 
in the SET and AI = e(F° - F^). The weights of the 
peaks are given by the initial weights of the eigenstates 



of Ha 



and \b\ . Measuring the charge m thus con- 



stitutes a perfect quantum measurement (Shnirman & 
Schon 1998). 

Next we analyze the effect of the perturbation E^^^pz in 
the master equation (|5.9|). It appears in both the coher- 



ent (LHS) and incoherent (RHS) parts of Eq. (5.9) and 
leads to the mixing. As usual the perturbation has the 
strongest effect when it lifts a (near) degeneracy, since 
in this case the eigenvectors within the degenerate sub- 
space may change substantially. Therefore we first treat 
the two near degenerate modes, i.e., we restrict ourselves 
to the two-dimensional subspace spanned by V!^'^{k = 0) 
and Vl^(k — 0), and check how the degeneracy between 
these two modes is lifted in second order of perturba- 
tive expansion. To account for this we approximate the 
diagonal part of the density matrix as 

/5diag(fc,t) ~ A°(fc,t)Fr(fc = 0) + A'{k,t)Vl\k - 0) , 

(5.26) 

(the LHS should be understood as an eight-column vector 
consisting of all matrix elements of Pdiag) for which we 
obtain an effective (reduced) master equation: 



dt [A^k) ) 



^^-''"'\A\k) ) 



(5.27) 



where 
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A"(fc) 



2t„ 



(5.28) 



The second term in ( |3.28 ) results from the perturbative 
expansion and, indeed, Ufts the degeneracy between the 
two modes. The mixing rate is obtained in second order 
perturbation expansion as 



-^int 



(5.29) 



In the approximation ( 5.26 ) we have A'{k,t) = 
J2n PnW- At t = these are occupation probabil- 
ities of the eigenstates of Hav, A°{k,t = 0) = |ap, 
A^{k,t = 0) = f or all k. It is straightforward to 
diagonalize Eq. ( 5.28 ) to find the eigenvalues. While 
the sum of the occupation probabilities is conserved, 
A°{k = 0,t) + A^{k = Q,t) = 1, the difference decays, 
A°lk = 0,t) - A^{k = 0,t) cx cxp(-</T„iix). Thus 
both occupations tend to 1/2 in the long-time limit. 
This implies that the evolution of P(rn, t) is given by 
Eq. (5.23) only for times t <C r^i-x, and in order to per- 
form a measurement one must have Tmix 3> T,„oas- As 



cx cos^ ?7av and t„ 



' mcas • 

2 • 2 

(X sm ' 



int ^ '/av 'mix -^int °^ Sm~77av, " 

good measurement may always be achieved by choosing 
tan^ ?7av small enough. 

Another effect of the mixing perturbation is that the 
dephasing rate, i.e., the rate of vanishing of the off- 
diagonal elements of the density matrix is changed as 



T^^ COS^ ?7av + l/2Tn 



(5.30) 



This is analogous to the relation between the dephas- 
ing ( 4.12D and relaxation (4.11) rates in the spin-boson 
model. It reflects the fact that the relaxation of the di- 
agonal elements of the density matrix leads also to ad- 
ditional suppression of the off-diagonal ones. Thus even 
when ?7av = 7''/2, the off-diagonal elements vanish with 
the rate l/2T,iiix. 

The long-time behavior of the qubit and detector, ex- 
cluding the period of the initial dephasing, is dominated 
by the two slowly decaying modes. In this regime we 
obtain from Eqs. (5.27,5.28) the reduced time evolution 
operator 



?7rcd(fc,t) = exp[M,ed(fc)i] ■ 



(5.31) 



Its Fourier transform Urcd{'ni,t) yields the expression for 
the distribution 



P(m,i) = (l, l)-C/rcd(m,i) 



\b\' 



(5.32) 



This Fourier transform can be performed analytically 
(Makhlin et al. 2000d), and we arrive at 



P(m,t)= ^ P{m-m\t) 



a'V2/rt 



v/27m 



(5.33) 



where 



m — Tt 



t 



+e-'"^— [\a\'S (m - Th) + |6|"J (m - Th)] 

and 



(5.34) 



(5.35) 



We observe that the solution is constructed from two 
delta-peaks, smeared by the convolution with the shot- 
noise Gaussian, and a plateau between them, which for 
|a;| < 1 is given by 



1 



AOL 'Yaiyi 

+ [1 + .x(|a|2 - l^n] h [t^JI-x^^ /v/T^} , (5.36) 

while Ppi = for |a;| > 1. Here /o, Ii are the modified 
Bessel functions. We work in the limit of weak qubit- 
detector coupling, where the Fano factors ( 5.2Cl| ) are close, 
and we denote them simply by /. For short times, t <^ 
Tjuix, the peaks dominate, while the plateau is low. At 
longer times the initial peaks disappear while the plateau 
is transformed into a single central peak (see Fig. p^ ). 

We complement the results for the charge distribution 
function P(m, t) by a derivation of the distribution fimc- 
tion for the values of the current. The measured quantity 
is actually the current averaged over a certain time in- 
terval At, i.e., /At = [m{t + At) - m{t)]/At = Am/ At. 
Accordingly the quantity of interest P{lAt,t) can be ex- 
pressed by the joint probability that m electrons have 
tunneled at time t and m -f Am electrons at a later time 
t + At: 

P{lAt,t) = ^ P(m, t:m + Am, t + At) 

711 

= Tt' [U{Am, At) U{m, i)/5o] • (5.37) 

m 

Here we made use of the Markovian approximation (see 
Section |V.El| ). The trace is taken over all degrees of free- 
dom except m, and also po refers to their initial state. 
In the two-mode approximation (5.27, 5T2^ ), sufficient at 
long times, we replace U again by C/rod and obtain 



P{lAut)^{\ l)-C/„d(Am,At) 
1 / p ^ |«P-|6|% ,,/. 



(5.38) 
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FIG. 19. Probability distribution P{lAt,t) of the current av- 
eraged over various time intervals At. The time axis is plotted 
on a logarithmic scale. 

The behavior of P(/At, is displayed in Fig. ^for var- 
ious values of At: (a) If the current is averaged over very 
short intervals, At <C r„ieas [Fig. ^ (a)], the detector 
does not have enough time to extract the signal from the 
shot-noise governed background, (b) An effective quan- 
tum measurement is achieved if Tmoas < At < Tnux- In 
this case the qubit-sensitive signal can be seen on top of 
the shot noise. Hence, as seen in Fig. pj| (b), the mea- 



sured value of Ja* is either close to eF" or to eF^ (5.19). 
The corresponding probabilities (the weights of the two 
peaks) are initially |ap and They change at longer 
times, t > Tinix, to a 1/2-1/2-distribution due to the mix- 
ing. A typical current pattern is a telegraph signal jump- 
ing between eF° and eF^ on a time scale Tj^ix- (c,d) If 
the current is averaged over longer times At > Tmix, 
the meter-induced mixing erases the information about 
the qubit's state. Several telegraph jumps can occur on 
this time scale, and for At 3> Tmix one measures only 
the time-averaged current between eF'' and eF^, com- 
pletely insensitive to the initial values. This is shown 
in Fig. |l^ (d), while Fig. |l^ (c) displays the crossover 
between (b) and (d). 

Notice, that since P(/At,t = 0) = P{m = IAt,t = At), 
the zero-time limits of the surfaces in Figs. ^ (a)-(d) 
are given by the charge distribution function plotted in 
Fig. [l|. 

We point out that in the Hamiltonian dominated 
regime the current in the SET is sensitive to the occupa- 
tion probabilities of the eigenstates of Tiav , rather than 
those in the basis of charge states. This may appear sur- 
prising, since the SET couples to the charge operator of 
the qubit. More precisely, the current is only sensitive to 



the expectation value of the charge in each eigenstate. As 
a consequence, at the degeneracy point Tyav = 7r/2, where 
the two eigenstates have the same average charge, both 
eigenstates lead to the same current in the SET and no 
measurement is possible. The measurement is effective 
only when charges in the eigenstates differ, tanryav ^ 1. 



F. Detector-dominated regime 



When the back-action dcphasing dominates over the 
average Hamiltonian, t^q 3> Ai?av, a perturbative anal- 
ysis in the charge basis is appropriate. In this basis the 
perturbation is the Josephson term of the qubit's Hamil- 
tonian Tictrb i-G., -^{i/'2)Bx<7x- Starting from the zeroth 
order, ^ Ej = 0, we obtain equations similar to 
Eqs. ( 5.18 , |3.21 ) with the replacements Ai?av B^^av 



and Er^^ Ei^t, where B^^av = B^ - 2{N)Eint is the 
charging energy of the "average" Hamiltonian. The anal- 
ysis of the diagonal and the off-diagonal modes is per- 
formed similarly. We get for the dephasing rate = 



. For the measurement time we reproduce Eq. (5.24). 
The dynamics of the two long-living diagonal modes can 
again be reduced to Eqs. ( p| ), 1 ^^ , (|^ and ( |5^ ) 
with 



EjT^o 



(5.39) 



This result is standard for the Zcno regime, i.e., the 
regime when coherent oscillations are overdamped by de- 
phasing (cf. Eq. 4.18). 

The condition for the Zeno regime given above re- 
quires a rather strong dephasing, such that t~q ex- 
ceeds both components of the qubit's Hamiltonian (^.3|), 
T^o ^ Ej,Bz,av The second part of this condition can 
be satisfied by tuning the qubit close to the degeneracy 
point. In contrast, in the Hamiltonian-dominated regime 
it is desirable for a good measurement to switch the qubit 
away from the degeneracy point. 

The long-time behavior of the charge and current 
distri butions, P( m,t) and P(/At,t), is again given by 
Eqs. (Ul imi) but with proper redefinitions of F , 



F , r~j^ etc. which now refer to the detector-dominated 
regime. In particular, the measurement provides infor- 
mation about the initial occupations of the qubit charge 
states rather than the eigenstates. 

To summarize we present a table with the main results 
for the two regimes: 
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regime 
leading term 


Hamiltonian 
dominated 
(coherent) 


fluctuation 
dominated 
(Zeno) 

SN{t)Snint 


largest energy 
preferred basis 
perturbation 


A^av > 

eigenbasis 

E\-nt sin ?7av Pr 


charge basis 

Ej <Tx 




'''ipO '^^^ '^av + 2 ''"mix 




_ 1 

r 

mcas 


7-i2 2 
OC E-^^ COS 77av 


7-l2 


mix 


4r(i;i„t sin?7av)^ 







Therefore only statistical (weak repeated) measurements 
are possible in this regime. 



G. Flux measurements 

To measure the quantum state of a flux qubit (for sim- 
plicity we consider here a single junction rf-SQUID), one 
needs a dissipative device sensitive to the magnetic flux. 
Natural candidates are dc-SQUIDs, i.e., superconduct- 
ing loops interrupted by two Josephson junctions. Cou- 
pling the qubit and the SQUID inductively, as shown in 



Fig. 20, one transfers part of the qubit 's flux into the loop 
of the SQUID and, thus, makes the SQUID sensitive to 
the state of the qubit. 

SQUIDs have been used as ultra-sensitive magnetome- 
ters for many years, and extensive literature covers the 
physics of the SQUIDs in much detail (see e.g. Tinkham 
1996, Likharev 1996). Here we mention only some of 
the recent activities and ideas which arose in connection 
with flux qubits. Two main strategies were proposed. 
The first one is to use underdamped dc-SQUIDs in the 
hysteretic regime (Mooij et al. 1999). This regime is re- 
alized when the SQUID is unshunted or the shunt re- 
sistance is large, Rs ^ {Ej/Ec) Rk, where Ej and Eq 
are the characteristic Josephson and charging energies 
of the SQUID. The higher the shunt resistance Rs, the 
lower is the noise of the SQUID in the superconducting 
regime (when no measurement is performed). Therefore, 
in spite of the permanent presence of the SQUID, the 
coherent dynamics of the qubit suffers only weak dephas- 
ing. To read out the state of the qubit, the current in the 
SQUID is ramped, and the value of current is recorded 
where the SQUID switches to the dissipative regime. As 
this switching current depends on the flux through the 
SQUID, information about the state of the qubit is ob- 
tained. Unfortunately, the current switching is a random 
process, fluctuating even when the external flux in the 
SQUID is fixed. For currently available system parame- 
ters this spread is larger than the difference in switching 
currents corresponding to the two states of the qubit. 



I* 




FIG. 20. Measurement setup for a flux qubit. The qubit (the 
rf-SQUID on the left side) is inductively coupled to the meter 
(the shunted dc-SQUID on the right side). 



The second strategy is to use overdamped SQUIDs, 
with Rs < {Ej / Ec) Rk- When the bias current exceeds 
the critical value, the voltage which develops across the 
shunt resistor depends on the external flux in the SQUID. 
Thus by measuring this voltage one learns about the state 
of the qubit. In this regime the principle of the mea- 
surement is identical to the one presented above for the 
SET. Recently Averin (20006) analyzed continuous (sta- 
tionary) measurements in this regime and obtained the 
input and output noise characteristics (cf. Section V.H) 
which determine the relevant time scales t^, t. 



mcasi 'mix- 



The main disadvantage of this strategy is that the SQUID 
induces dephasing during the periods of coherent manip- 
ulations when no measurement is performed. The ques- 
tion still remains to be settled whether a reasonable com- 
promise between the underdamped and the overdamped 
limits can be found. 



H. Efficiency of the measuring device 

Recently, several devices performing quantum mea- 
surements have been analyzed. Apart from SETs in 
the sequential tunneling regime (Shnirman & Schon 
1998, Makhlin et al. 2000 d, Korotkov 2000 &, Devoret & 
Schoelkopf 2000) it includes SETs in the co-tunneling 
regime (Averin 2000 a, Maassen van den Brink 2000), su- 
perconducting SETs (SSETs) and dc-SQUIDS (Averin 
2000 &), as well as quantum point contacts (QPC) (Ko- 
rotkov & Averin 2000). All these devices are based on 
the same basic idea: they are dissipative systems whose 
response (conductance, resistance) depends on the state 
of a qubit coupled to them. 

The efficiency of a quantum detector has several as- 
pects. From a practical point of view the most important 
is the ability to perform a strong, single-shot measure- 
ment which requires that the mixing is slower than the 
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read-out, Tmix ^ Tmeas- Another desired property is a 
low back-action noise of the meter in the off-state, which 
can be characterized by a corresponding dephasing rate. 

A further important figure of merit is the ratio of the 
dephasing and measurement times. Quantum mechan- 
ics demands that a quantum measurement should com- 
pletely dephase a quantum state, i.e. < r^cas- For 
the example of a SET coupled to a charge qubit the de- 
phasing time ( 5.22| ) is smaller (or even much smaller) 
than the measurement time ( |5.24| ). This means that 
the information becomes available later than it would 
be possible in principle. In this sense the efficiency of 
the SET in the sequential tunneling regime is less then 
100% (Shnirman & Schon 1998, Korotkov 1999). The 
reason for the delay is an entanglement of the qubit 
with microscopic degrees of freedom in the SET. To il- 
lustrate this point consider a situation where the initial 
state of the system (a|0) -I- 6|1)) \x) \m = 0) evolves into 
o|0) Ixo) \iTio) + Ixi) where |x) stands for the 

quantum state of the uncontrolled environment. One 
can imagine a situation when mg = mi, but \xo) and 
Ixi) are orthogonal. Then the dephasing has occurred 
but no measurement has been performed. It is interest- 
ing to note that the ratio r^/rmcas grows if the SET is 
biased in an asymmetric way, creating a strong asym- 
metry in the tunneling rates, e.g. Fl <C Fr. In other 
measurement devices the ratio r^p/Tmcus may be close 
to 1. This includes quantum point contacts (Korotkov 
1999, Korotkov & Averin 2000) as well as SETs in the co- 
tunneling regime (Averin 2000 a, Maassen van den Brink 
2000). The common feature of these three examples is 
that the device consists of one junction or effectively re- 
duces to it. It should be kept in mind, however, that 
a large ratio of dephasing and measurement times is not 
the only figure of merit. For instance, in a SET in the co- 
tunneling regime the current is low and more difficult to 
detect as compared to a SET in the sequential tunneling 
regime. 

The ratio of dephasing and measurement times 
Ti^/Tmoas has been analyzed also in the framework of 
the theory of linear amplifiers (Averin 20006, Devoret 
& Schoelkopf 2000). It can be expressed in terms of the 
noise characteristics of the amplifier which, in turn, de- 
termine the sensitivity of the device. In this framework, 
one considers a detector with the output signal / and in- 
put signal 0, which is coupled to an observable Q of the 
detector via a term (f>Q. The input signal causes a vari- 
ation of the output, which can be characterized by the 
linear response coefficient A = d{I)/d(p. Note that usu- 
ally one operates in a dissipative, nonequilibrium regime. 
When used as a quantum detector coupled to a qubit the 
input variable is cf) (x a^, and the coupling is ccr^Q (with 
c being a coupling constant). We first consider a situ- 
ation where tunneling between the qubit 's basis states 
is suppressed, i.e., TYctri ~ —Bz(7z/2. In this case the 



fluctuations of Q dephase the qubit with rate 



(5.40) 



The symbols Sq and Si (introduced below) stand for 
the noise power of the corresponding observable if the 
amplifier is decoupled from the qubit: Sq = 2 (Q^) and 
Sj = 2 (/^). A white spectrum is assumed at the relevant 
frequencies. The two basis states of the qubit produce 
output signals /"^^ — I ± AI/2, differing by 



AI = 2cA 



(5.41) 



They can be distinguished after the measurement time 



_ (Air 

mcas 

Hence, the two times are related by 
_ SqSi _ SqS^ 



(5.42) 



(5.43) 



In the last form of Eq. ( 5.43| ) S^ = Si/X^ is the output 
noise in terms of the input, i.e., the noise which should 
be applied to the input to produce the noise Si at the 
output. 

If the tunneling is turned on, Tictri — — ■^A i?(cos r] (Tx + 
sin?7 az), both the measurement rate (5.42) and the rate 
of pure dephasing (5.40) acquire an additional factor 
cos'^T^, while their ratio ( 5.43| ) persists. Apart from that, 
the finite tunneling introduces mixing, with rate 



^S'q sin f] 



(5.44) 



For the particular case of the SET, which motivates the 
notations, / is the transport current, Q = eN the charge 
of its central island, and (f) the gate potential externally 
applied to the middle island of the SET (0 consists of 
Vat and a contribution of the qubit, see Eq. The 
coupling constant is c = i^int / e and the linear response is 
given by A = A//(2c) = e^ST/{2Eint) with 5T defined by 



Eq. ( 5.35| ). In this case, Eqs. ( ^.4C| , p3.42 ) are consistent 
with what has been described before in Eqs. (5.22, ^.25 ). 

The quantity y^SqS^ is proportional to the "noise en- 
ergy", discussed by Devoret & Schoelkopf (2000), mea- 
sured in units of the energy quanta at the given frequency. 
For the noise spectra Sq and S^ one can obtain an in- 
equality, similar to the Heisenberg uncertainty principle 
(Braginsky & KhaliU 1992): SqS^ > Ti^ . Indeed, by 
virtue of Eq. (5.43) this relation coincides with the con- 
straint 



"^if — "^mcas 



(5.45) 



32 



Moreover, one can show (Braginsky & Khalili 1992, 
Averin 2000?j)-jthat a stronger inequality holds. Namely, 
the quantity Il3 



4>Q) 



(5.46) 



is limited by e > h/2. Here = S'/q/A, where 

SiQ = 2 J dt'{I{t)Q{t')) characte rizes cross-correlations 
between the input and the output. Thus one has 



\h/2j + 



(5.47) 



The optimization of this ratio requires that the detec- 
tor reaches the quantum limit of sensitivity, e ~ ^/2, 
and the cross-correlations vanish, ReS'^g = 0. Averin 
(2000 &) made the observation that the quantum limit is 
(nearly) reached in several measurement devices: For a 
quantum point contact, an overdamped dc-SQUID, a re- 
sistively shunted superconducting SET, or a normal SET 
in the co-tunneling regime (Averin 2000 a) the following 
two relations hold under certain conditions: 



SrSr 



\S_ 



IQ\ 



and 



IQ 



(5.48) 



(5.49) 



They im mediately imply that e = fi/2. Note, that the 
relations ( 5.48 5.4E ) are requirements for a quantum lim- 
ited measurement device. They are not valid in gen- 
eral. E.g., the near equality in ( 5.48 ) should in general 
be replaced by the > sign. Even for the detectors where 
they were found they break down (making e > ?i/2), e.g., 
at finite temperatures or at bias voltages close to the 
Coulomb blockade threshold (Averin 2000 a). 

As for Eq. ( |3.4g| ) , by definition the response coefhcient 
is X = -ij dt'e{t - t'){\I{t),Q{t')]) = 2Im / dt'e{t ~ 
t'){I{t)Q{t')). Hence, Eq. ( ^.49D implies a vanish- 
ing reciprocal response coefficient A' = 2Im J dt'6{t — 
t'){Q{t)I{t')) . While in equilibrium Onsager's relations 
imply A = —A', such an asymmetry can arise in a 
nonequilibrium stationary state. It is a characteristic 
feature of linear measuring devices (Braginsky & Khalili 
1992). 

As for the cross-correlations, ReS'^g, they vanish, e.g., 
in a dc-SQUID with two identical Josephson junctions, 
a QPC, symmetrically coupled to a quantum dot, or a 
SSET (Averin 20006). 



^^This quantity is denoted as "energy sensitivity" by Averin 
(20006), but Devoret & Schoelkopf (2000) use this term for a 
different quantity. 



I. Statistics of the current and tlie noise spectrum 

In previous subsections we have discussed the statis- 
tics of the charge m, which passed through the detector, 
and the corresponding current / after a measurement has 
been started. In this subsection we complement this dis- 
cussion by investigating the noise properties of the charge 
and the current in the stationary regime, i.e., a long time 
after the qubit-detector coupling was turned on. The 
noise spectrum reflects the intrinsic properties of the sys- 
tem of qubit and meter and their coupling and depends 
on the corresponding time scales. We obtain additional 
evidence for the telegraph behavior of the current. 

We study first the charge-charge correlator, which is 
derived from the j oint p robability distribution of charges 
at different times (^.37|) , and obtain 



{m{ti) m{t2)) 



U{k,M)dkU\k,t^^ 



(5.50) 



k=0 



where = — and t = {ti+t2)/2. We employ again 
the two- mode approximation and use Eqs. (5.28,5.31) to 
arrive, after taking time derivatives, at 



^2 ^g-AtA 



+ (• 



(5.51) 



Here T^/^ = T ±5T/2 are the tunneling rates (^.35|) for 
the t wo q ubit states. The Fano factor / is defined after 
Eq. ( 5.34 ). In the stationary limit, i 3> Tmix, we thus 
obtain for the noise spectrum 



Si{uj) = 2eVf + 



1 -t- wV^. 



(5.52) 



The first, tj-independent term corresponds to the shot 
noise, while the second term originates from the tele- 
graph noise. The ratio of the telegraph and shot noise 
amplitudes at w = is 



telegraph 



shot 



2/f 



(5.53) 



Note that the telegraph noise becomes noticeable on top 
of the shot noise in the parameter regime of an effective 
quantum measurement (rmix 3> Tmoas)- 

While the shot-noise contribution reflects intrinsic 
properties of the detector (SET), the telegraph noise 
characterizes the qubit. This structure of the noise of 
the output signal is quite general (Korotkov & Averin 
2000, Averin 20006, Korotkov 2000 a). For the qubit cou- 
pled to the meter, the following relation for the noise of 
the output signal / was derived (using linear response 
and certain other approximations): 



Sl{uj) = S%UJ) + C^X^ SaA^) 



(5.54) 
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Here the notations introduced in Section V.H are used, 



with the only difference that Sj now denotes the noise 
in the presence of the qubit. The first term, Sj, in the 
RHS of Eq. (5.54) represents the noise of the meter de- 
coupled from the qubit (^shot for a SET). The second 
term arises from the qubit and is governed by the dy- 
namics of the its density matrix. In the preferred basis 
(cf. Section V.D-V.F) the diagonal elements decay to 



their stationary values on the mixing time scale, while 
the off-diagonal elements precess with frequency AE/h 
and decay within the dephasing time. Depending on the 
qubit 's Hamiltonian, both or one of these processes con- 
tribute to the dynamics of (7z{t). Hence, the qubit's con- 
tribution, S'o-^(w) = 2/ dt{az{0)az{t)) exp{—iujt), has a 
"coherent" peak at the qubit's eigenfrequency, AE/h, of 
width T~^, and a "telegraph" peak at zero frequency with 
width r^i^^. 

To derive the noise spectrum of the qubit under the 
influence of the detector, one can use a Bloch-type mas- 
ter equation for the joint density matrix [e.g., Eq.( |5.9| ) 
at fc = for a SET]. Under certain conditions, namely, if 
the noise spectrum Sq^lo) is white in an interval of fre- 
quencies at least up to the qubit's level spacing AE (in 
the SET it is white only up to Fr + Fl) a simpler set of 
equations is sufficient (Averin 2000 &): 



Pqo 



Poi = {iBz - To)poi - ^B^ipoo - pii) . 



(5.55) 



Here pij is the qubit's density matrix in the eigenbasis 



of. 



c^Sq/Ti^ describes the back-action ( ^.40 ) 



and Bx and Bz form the Hamiltonian of the qubit (p. 3). 
These equations were used, for instance, for the analysis 
of a quantum point contact coupled to a double quantum 
dot (Gurvitz 1997, Korotkov & Averin 200 0). 

At zero bias, Bz = 0, one can solve Eq. ( 5.55| ) exactly, 
to obtain (Korotkov & Averin 2000) 



(5.56) 



In the Hamiltonian-dominated limit, t^^q <C Bx, at this 
bias point no telegraph peak appears. Still, the qubit 



contributes to the noise via the last term in Eq. (5.54) 
This coherent peak at B^ has the height 



{Alf 



(5.57) 



Here AI is the difference in the output current ( ^.41 ) for 
the two eigenstates of CTz (charge states). Recalling that it 
is related to the measurement rate (5.42) in the detector- 
dominated regime, one concludes that the signal-to-noise 
ratio in Eq. ( 5.54| ) is limited, by virtue of Eq. ( 5.45| ), as 



5° 



= 4 



T<pO 



< 4 . 



(5.58) 



Thus, the requirements for the observation of coherent 
oscillations in the noise spectrum of the measuring de- 
vice is directly expressed by the ratio of dephasing and 
measurement time. 

In the opposite, detector-dominated regime B^ <SC 
T^Q , the noise spectrum ( 5 .561 ) exhibits a telegraph- noise 



Lorentzian ( 5.52| ) at low frequencies, 4Ti„ix/(w^r,^^ix + 1) 



where r„ 



Bi 



At a general bias point Bz ^ in the Hamiltonian- 
dominated regime the dynamics of <7z exhibits features 
at both frequencies, a; = and to = AE, and its noise 
spectrum has two peaks. This was shown numerically 
by Korotkov & Averin (2000) and it is consistent with 
the result (5.52). A higher-order analysis of the master 
equation (5.£), which accounts for the corrections to the 



oscillating modes (5.21) due to the mixing, allows us to 



obtain the coherent peak which is, however, suppressed 
due to the strong dephasing (cf. Eq. 5.58 ). 



J. Conditional master equation 

We finally would like to comment on recent studies 
(Korotkov 1999, Korotkov 2000a, Goan et al. 2000) where 
the so-called selective or conditional approach, popular in 
quantum optics, was employed. Its purpose is to develop 
a framework for the analysis of statistical properties of 
qubit-related quantities conditioned on the dynamics of 
the output signal, I{t), at earlier times. In other words, 
the goal is to predict the outcome of a measurement of 
the qubit's state at time t given the results of the cur- 
rent read-out at t' < t. A related problem is to produce 
typical, fluctuating current-time patterns I(t) which one 
can encounter in a given experiment. 

These problems can be addressed using the proper 
master equation for the coupled system (e.g., Eq. |3.9[ ). 
Monitoring the current amounts to repeated measure- 
ments of the charge m, at sufficiently short time- intervals 
At. Since a closed master equation can be formulated 
which involves only the m-diagonal entries of the density 
matrix p{m), the effect of a measurement with the result 
Too amounts to choosing the corresponding p(rno)Smmo 
as the new density matrix (properly rescaled to ensure 
the normalization). Thus a simulation of the evolution 
proceeds as follows: From the the master equation with 
initial matrix pSmo &t t — one obtains p{Am, At) at the 
time of the first read-out. Then, to simulate an experi- 
ment, one selects a "measured" value of Am with the cor- 
responding probability P(Ato, At) = tr p{Am, At) (cf. 
Section |V.E|) and uses the corresponding density matrix 
as initial value for the further evolution during the next 
time interval. Repeating this step many times one pro- 
duces a typical dependence of m{t) and a density ma- 
trix Pcond{t) which can be used to study the conditional 
statistics of further measurements. 
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This procedure can be simphfied if sufficiently frequent 
read-outs are performed. An expansion in At allows one 
to present the step-by-step evolution of the density ma- 
trix pcond as a continuous process (cf. Korotkov 20006): 



d 
di 



p = -i[HctrhP] - ilU^z, [(^z,p]] 
AI 



1 

4 "fi^ 



+ ^m-I][{a,„p}-2tr{a,J)p] , (5.59) 



where AI is defined in Eq. (5.41). The measured value 
of the output signal, I{t) — Am/ At, as described 
above, should be chosen randomly with the distribution 
P{A'm, At). Using the properties of P{m,At) at short 
At, one can further simplify this procedure and choose 
I(t) as 



I{t)=I+^tr{a,J) + SI{t) , 



(5.60) 



where 61{t) is a random quantity with the noise prop- 
erties identical to those of the current in the detector 
decoupled from the qubit. 



We can menti on th at, similar to Eq. ( 5.55 ), in the 
derivation of Eq. (5.59) we assumed the back-action noise 



of the detector to be white. Furthermore, 



is the de- 



phasing rate due to the environment and accounts for 
non-ideality of the detector (see Section |VJ|. After av- 
eraging over the detector dynamics, I{t), Eq. ( 5.59 ) is 
recovered. 



The equations (5.5£, 5.60) form a Langevin-type evo- 
lution equation which can be used to produce typical 
outcomes of the measurement and to study statistics of 
the qubit conditioned on these outcomes. 



VI. CONCLUSIONS 

Josephson junction systems in a suitable parameter 
range can be manipulated in a quantum coherent fashion. 
They are promising physical realizations (the hardware) 
of future devices to be used for quantum state engineer- 
ing. We discussed their modes of operation in different 
designs (in the charge and the flux dominated regimes), 
the constraints on the parameters, various dephasing ef- 
fects, and also the physical realization of the quantum 
mechanical measurement process. We pointed out the 
advantages of these nano-electronic devices as compared 
to other physical realizations. 

We add a few remarks and comparisons. First, in order 
to demonstrate that the constraints on the circuit param- 
eters, which were derived in previous sections, can be met 
by available technologies, we summarize them here and 
suggest a suitable set. 

(i) Necessary conditions for a Josephson charge qubit are: 
A > Eq ^ Ej , k-oT. The superconducting energy gap 



A has to be chosen large to suppress quasiparticle tun- 
neling. The temperature has to be low to assure the 
initial thermalization, k^T <C Eq, TilolCi a-nd to reduce 
dephasing effects. A sufficient choice is fceT ~ Ej/2, 
since further cooling does not reduce the dephasing (re- 
laxation) rate in a qualitative way. (Of course, it does 
so far from the degeneracy point, i.e., for rj = 0, or if we 
switch off the Hamiltonian, Tictri = 0. However during 
manipulations Ej is the typical energy difference and sets 
the time scale for both the manipulation times and the 
dephasing.) 

As an explicit example we suggest the following pa- 
rameters (the circuit parameters of Nakamura et al. 
(1999) are in this regime) and estimate the correspond- 
ing time scales: We choose junctions with capacitance 
Cj = 10~^^ F, corresponding to a charging energy (in 
temperature units) Eq ^ 1 K, and a smaller gate ca- 
pacitance Cg — 0.5 • 10^^^ F to reduce the couphng to 
the environment. Thus at the working temperature of 
T = 50 mK the initial thermalization is assured. The 
superconducting gap has to be slightly higher, A > Ec- 
Thus aluminum is a suitable material. We further choose 
Ej = 100 mK, i.e., the time scale of one-qubit operations 
is rip^ = h/E] - 10-1" s. 

(ii) A realistic value of the resistor in the gate voltages 
circuit is i? ~ 50f2. Its voltage fluctuations limit the de- 



phasing time ( |4.12| ) to values of order ~ 10 s, thus 

allowing for jzlj/top'' ^ 10^ coherent manipulations of a 
single qubit Ej. 

(iii) To assure sufficiently fast 2-bit operations we choose 
for the design of Fig. I L - 10 nH and Cl w Cj. Then 
the 2-bit operations are about 10^ times slower than the 
1-bit operations and, accordingly, their maximum num- 
ber is reduced. 

(iv) The quantum measurement process introduces addi- 
tional constraints on the parameters, which can be met 
in realistic devices as demonstrated by the following con- 
cise example. The parameters of the qubit are those 
mentioned before. For the junction and gate capaci- 
tances of the normal tunnel junctions of the SET we chose 
Ct = 1.5-10-1^ F and Cf^"^ = 0.5 -lO"!^ F, respectively, 
and for the coupling capacitance between SET and qubit: 
Cint = 0.5 • 10-1^ F. We thus obtain: Ec,set ~ 25 K, 
-E'c.qb ~ 1 K, Eint ~ 0.25 K (for precise definitions see 
Appendix |C|) . We further take ng,qb — 0.25, iVg.sEx = 
0.2 and pL = -PR = eVir/2 = 24 K and = or = 0.03. 
This gives B;, « 2 K and Fl = 1.8 K and Fr = 7.8 K 



^^This may be overly optimistic, and indicating that other 
sources of dephasing need to be considered as well. For in- 
stance, at these slow time scales the background charge fluc- 
tuations may dominate. We also note that in the experiment 
of Nakamura et al. (1999) a stray capacitance in the probe 
circuit, larger than Cg, renders the dephasing time shorter. 
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(note that due to the gate voltage the apphed transport 
voltage is split asymmetrically). Thus the measurement 
time in this regime is Tmcas ~ 1-5 ■ 10^^ 
choice of parameters 



1/30 



s. Since for this 
^ ~ 4.0 • 10-3 K < 



the "Hamiltonian dominated" regime is realized. Assum- 
ing E,^ = 0.1 K we obtain r„iix « 0.7 • 10"'^ s. Thus 
Tmix/Tmcas ~ 45 and the separation of peaks should be 
observable early enough before the mixing dominates. In- 



deed, the numerical simulation of the system (5.9) for 
these parameters shows almost ideal separation of peaks 
(see Fig. |l5|). On the other hand, for Ej = 0.25 K we 
obtain T,„ix/Tnioas ~ 7. This is a marginal situation. The 
numerical simulation in this case (see Fig. |l^) shows that 
the peaks first start to separate, but, soon, the valley be- 
tween the peaks fills in due to the mixing transitions. 

The requirements on the parameters of a flux qubit 
circuit can be summarized in a similar way. First, the 
parameters should be chosen to allow the reduction of 
the double-well potential to the two ground states form- 
ing a two-state quantum system. It is also desirable that 
these two basis states have macroscopically different flux 
or phase configurations. The double well is formed by 
joining either several Josephson junctions (Mooij et al. 
1999) or a Josephson junction and an inductive term of 
similar strength (i?j ^ ^l/An^L) in an rf-SQUID (Fried- 
man et al. 2000)). Since the level spacing within each 
well is of order Eq ~ ^/EJEq and the barrier height of 
order Ej, all these requirements can be satisfied by 'clas- 
sical' Josephson junctions with Ej ^ Ec- Furthermore, 
the asymmetry of the double well Bz, which is controlled 
by external fluxes, and the tunnel splitting should be 
smaller than Eq in order to suppress leakage to higher 
states. Finally, the temperature should be low enough to 
allow the initialization of the qubit 's state and to ensure 
slow dephasing. In summary, the following conditions 
have to be satisfied: k^T < B^, ^/EjEc < Ej. As 
pointed out above, it is sufficient to choose k-^T ^ B^jl- 
These requirements can be satisfied, e.g., by the param- 
eters of the rf-SQUID used by Friedman et al. (2000) 
[Ei « 70 K, ~ 1 K, B^ w 0.1 K, and T w 40 mK) 
or by similar values discussed by Mooij et al. (1999) 
{Ei - 10 K, Ec - 0.1 K, B^ - 50 mK, and T - 30 mK). 

At present, the most advanced quantum manipulations 
of a solid state system, i.e., the coherent oscillations ob- 
served by Nakamura et al. (1999), have been demon- 
strated for a Josephson charge qubit. But flux systems 
might soon catch up. In the long run, it is not clear 
whether charge or flux systems will bring faster progress 
and further reaching demonstrations of complex quantum 
physics. In fact, a combination of both appears feasible 
as well. Therefore we compare shortly the properties of 
the simplest charge and flux qubits. 



A very important quantity is the phase coherence time 
Tip, which has to be compared to the typical operation 
time scale of the qubit's dynamics Top. While effects of 
various dissipative mechanisms have been estimated the- 
oretically, further experimental work is needed to under- 
stand dephasing in charge and ffux qubits. One poten- 
tially dangerous source of dephasing, the coupling to the 
external circuit, can be described by an Ohmic oscillator 
bath. If this contribution to the dephasing dominates, 
the above-mentioned ratio of times is determined by the 
dimensionless parameter Top/r^ ~ a (see Section IV). 
For "unscreened" charge qubits it is of order of a w 10^^, 
but it can be substantially reduced by the 'screening ra- 
tio' of capacitances (see Section LV). Putting in numbers 
corresponding to those of Nakamura et al.'s (2000) ex- 
periment (without probe junctions), we estimate to 
reach several tens of microseconds. Coherent oscillations 
for about 5 ns have been observed already in this system, 
in spite of the presence of a non-ideal detector limiting 
the phase coherent evolution by quasiparticle tunneling 
in the probe junction and also providing a strong cou- 
pling to the external circuit. The phase coherence time 
should be compared to the qubit operation time scale of 
Top w 10 to 100 ps. 

For flux qubits considered by Mooij et al. (1999) the 
current circuit may produce a. ^ 10~^ and the relaxation 
times of order 10 /iS, while other dephasing mechanisms 
studied would destroy coherence after times of order of 
hundreds of microseconds or longer. At the same time, 
the qubit level spacing sets the fastest operation time to 
Top « 0.1 to 1 ns. 

A major source of errors and dephasing for all charge- 
degrees of freedom are the fluctuations of the 'offset 
charges'. They arise due to charge transfers in the sub- 
strate, e.g., between impurity sites, and are detrimental 
for many of the potential applications of single-charge 
systems. Fortunately, they occur typically on long time 
scales, and may not take place during a 'computation', 
i.e., a series of coherent manipulations. Similarly for 
flux qubits nuclear spins provide random magnetic fields. 
Also these fields change only on a long spin-relaxation 
time scale and cause no dephasing in shorter 'computa- 
tions'. 

Another important point is the efficiency of quan- 
tum detectors used to read out the charge or flux 
state. Estimates show that the newly developed rf-SETs 
(Schoelkopf et al. 1998) should make single-shot charge 
measurements possible in principle {Tmix ^ Tmcas)- On 
the other hand, the flux read-out with a SQUID is far 
from this goal and averaging over a large number of mea- 
surements is needed (van der Wal et al. 2000). 

As for the experimental achievements, in both charge 
and flux systems superpositions of basis (charge or flux) 
states have been seen, and the validity of the two-state 
model has been conflrmed spectroscopically (Bouchiat 
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1997, Nakamura et al. 1997, Friedman et al. 2000, van der 
Wal et al. 2000). As for the coherent oscillatfons between 
the basis states, so far they were seen only in the charge 
design (Nakamura et al. 1999). In Appendix ^ we fur- 
ther compare Josephson devices to some other systems 
suggested as physical realizations of qubits. 

To conclude, the fabrication and controlled coherent 
manipulations of Josephson junction qubits are possible 
using present-day technologies. In these systems fun- 
damental properties of macroscopic quantum-mechanical 
systems can be explored. First experiments on elemen- 
tary systems have been performed successfully. More 
elaborate designs as well as further progress of nano- 
technology will provide longer coherence times and allow 
sequences of coherent manipulations as well as scaling to 
larger numbers of qubits. The application of Josephson 
junction systems as elements of a quantum computer, 
i.e., with a very large number of manipulations and large 
number of qubits, will remain a challenging issue. On 
the other hand, many aspects of quantum information 
processing can initially be tested on simple circuits as 
proposed here. We also expect further spin-offs, once 
the techniques of quantum state engineering get further 
developed. 

We have further shown that a dissipative quantum sys- 
tem coupled to a qubit may serve as a quantum measur- 
ing device in an accessible range of parameters. Explic- 
itly we studied a single-electron transistor coupled capac- 
itively to a charge qubit. We have described the process 
of measurement by deriving the time evolution of the re- 
duced density matrix of the coupled system. We found 
that the dephasing time is shorter than the measurement 
time, and we have estimated the mixing time, i. e. the 
time scale on which the transitions induced by the mea- 
surement occur. Similar scenarios are discussed for flux 
qubits measured by a SQUID coupled to it inductively. 
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APPENDIX A: AN IDEAL MODEL 

The model Hamiltonian 

Quantum state engineering requires coherent manipu- 
lations of suitable quantum systems. The needed quan- 



tum manipulations can be performed if we have suffi- 
cient control over the fields and interaction terms in the 
Hamiltonian. As an introduction and in order to clar- 
ify the goal, we present here an ideal model Hamiltonian 
and show how the necessary unitary transformations can 
be performed. We can add that the Josephson junctions 
devices discussed in this review come rather close to this 
ideal model. 

As has been stressed by DiVincenzo (1997, 2000) any 
physical system which is considered as a candidate for 
quantum computation, and similar for alternative appli- 
cations of quantum state engineering, should satisfy the 
following criteria: (i) First of all, one needs well-defined 
two-state quantum systems (or quantum systems with a 
small number of states) . This implies that higher states, 
present in most real systems, must not be excited during 
manipulations, (ii) One should be able to prepare the 
initial state of the qubits with sufScient accuracy, (iii) A 
long phase coherence time is needed, sufficient to allow 
for a large number (e.g., > 10'*) of coherent manipula- 
tions, (iv) Sufficient control over the qubit 's Hamiltonian 
is required to perform the necessary unitary transfor- 
mations, i.e., single-qubit and two-qubit logic operations 
(gates). For this purpose one should be able to control 
the fields at the sites of each qubit separately, and to cou- 
ple qubits together in a controlled way, ideally with the 
possibility to switch the inter-qubit interactions on and 
off. In physics terms the two types of operations allow 
creating arbitrary superpositions and non-trivial coupled 
states, such as entangled states, respectively, (v) Finally, 
a quantum measurement is needed to read out the quan- 
tum information, either at the final stage or during the 
computation, e.g., for the purposes of error-correction. 

We consider two-state quantum systems (e.g./ spins), 
or systems which under certain condition effectively re- 
duce to two-state systems (charge in a box or flux in a 
SQUID near degeneracy points). Any single two-state 
quantum system can be represented as a spin-1/2 parti- 
cle, and its Hamiltonian be written as 7i(t) = —^B{t)&. 
Here a^.y.z are Pauli matrices in the space of states 
It) — (p) and II) = (°), which form the basis states 
of a physical quantity (spin, charge, flux, . . . ) which is 
to be manipulated. The "spin" is coupled to an effective 
"magnetic fleld" B. In an alternative notation used for 
two-state quantum systems the components Bz and B^.y 
correspond to an external bias and a tunneling ampli- 
tude, respectively. Full control of the quantum dynam- 
ics of the spin is possible if the magnetic field B{t) can 
be switched arbitrarily. In fact, arbitrary single-qubit 
operations can be performed already if two of the field 
components can be controlled, e.g., 

Ktri(i) = ~^Bz{t)az - \B,{t)(T, . (Al) 

If all three components of the magnetic field can be con- 
trolled the topological (Berry) phase of the system can 
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be manipulated as well (Falci et al. 2000). 

In order to manipulate a many-qubit system, e.g., to 
perform quantum computing, the magnetic field at the 
site of each spin has to be controlled separately. In ad- 
dition, one needs two-qubit (unitary) operations which 
requires controlling the coupling energies between the 
qubits. For instance, a system with the following model 
Hamiltonian would be suitable: 



1 ^ 



(A2) 



where a summation over spin indices a,b = x,y, z is im- 
plied. In ( [A^ ) a general form of coupling is presented, 
but simpler forms, such as pure Ising zz-coupling, XY-, 
or Heisenberg coupling are sufficient. 

The measurement device, when turned on, and residual 
interactions with the environment are accounted for by 
extra terms 7imeas(i) and Tires, respectively: 



(A3) 



During the manipulations the meter should be kept in the 
off-state, Timcas = 0. The residual interaction Tires leads 
to dephasing and relaxation processes. It has to be weak 
in order to allow for a series of coherent manipulations. 

A typical experiment involves preparation of an ini- 
tial quantum state, switching the 'fields' B{t) and the 
coupling energies J^l{t) to effect a specified unitary evo- 
lution of the wave function, and the measurement of the 
final state. 



Preparation of the initial state 

The initial state can be prepared by keeping the system 
at low temperatures so that it relaxes to the ground state. 
This is achieved by turning on a large value of ^ k^T 
for a sufficiently long time while Bx{t) — By{t) = 0. 
Then the residual interaction. Tiros relaxes each qubit 
to its ground state, |t). Switching B2{t) back to zero 
leaves the system in a well-defined pure quantum state. 
If Ti = 0, there is no further time evolution. 



Single-qubit operations 

A single-bit operation on a given qubit can be per- 
formed, e.g., by turning on Bx{t) for a time span r. As 
a result of this operation the quantum state evolves ac- 
cording to the unitary transformation 



Ux {a) = exp 



iBxT&x 

2h 



(A4) 



where a — Bxt/Ti. Depending on the time span r, an 
a = TT- or an a = 7r/2-rotation is performed, producing 



a spin flip (NOT-operation) or an equal-weight superpo- 
sition of spin states. Switching on Bz{t) for some time 
produces another needed single-bit operation: a phase 
shift between |t) and ||): 



(/3) = exp 



2h 
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(A5) 



where (3 = B^r/h. With a sequence of these x- and z- 
rotations any unitary transformation of the qubit state 
(single-qubit operation) can be performed. There is no 
need to turn on B,,. 



Two-qubit operations 

A two-bit operation on qubits i and j is induced by 
turning on the corresponding coupling J*-' (t) . For in- 
stance, for the XF-coupfing, J^f^a-aO-i, — ■P^{ax^x + 
ayCTy), in the basis Ititj), IT4i), liitj), \Ulj) the result 
is described by the unitary operator 



I I Q\ 

cos 7 isin7 

jsin7 cos 7 

V 1 / 



(A6) 



with 7 = 2J^^T/h. For 7 = tt/2 the operation leads to 
a swap of the states |tiij) and |iitj) (and multiplication 
by i), while for 7 = 7r/4 it transforms the state |tiij) 
into an entangled state ^ (ITiij) + * litti))- 

We note that apart from the sudden switching of 
Bl^^{t), J^^t), discussed above for illustration, one can 
also use other techniques to implement single-bit or two- 
bit operations. For instance, one can induce Rabi oscil- 
lations between different states of a qubit or a qubit pair 
by ac resonance signals; or perform adiabatic manipu- 
lations of the qubits' Hamiltonian to exchange different 
eigenstates (with occupations remaining unchanged) . We 
discuss some of these methods for particular physical sys- 
tems in Sections || and III . 



APPENDIX B: QUANTUM LOGIC GATES AND QUANTUM 
ALGORITHMS 

In Appendix ^ and Sections ||, |ll| we showed how el- 
ementary quantum logic gates can be realized by simple 
manipulations of concrete physical systems. Details such 
as the application of a magnetic field pulse or the type of 
the two-qubit coupling depend on the specific model. On 
the other hand, quantum information theory discusses 
quantum computation in realization-independent terms. 
For instance, it is customary to build quantum algorithms 
out of specific, 'standard' single- and two-qubit gates, 
some of which will be discussed below. Hence, one needs 
to know how to express these standard gates in terms 
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of the elementary operations specific for a given physical 
model. Furthermore, one may be interested in an op- 
timized implementation, in terms of time, complexity of 
manipulations, or amount of additional dissipation. Here 
we give several examples of standard gates and quantum 
algorithms and cover optimization issues later. 



Single- and two-qubit gates 

The quantum generalization of the NOT gate, 

NOT : 



1 

1 



(Bl) 



permutes the basis vectors |0) |1) and |1) |0). It can 



be performed as the x-rotation (A4) with the time span 
corresponding to a = tt (up to an unimportant overall 
phase factor): Ux{ct = tt) = i ■ NOT. Unlike to classical 
computation, in quantum logic there exists a logic gate, 
called vNOT, which when applied twice produces the 
NOT-gate: 



VNOT = i 



(B2) 



This gate0 is obtained by an x-rotation (A4) with a = 
tt/2, more precisely Ux{a = n/2) = i^/^ • VNOT. 

Another important, essentially quantum mechanical 
single-bit operation is the Hadamard gate: 



H 



V2\l -I 



(B3) 



It transforms basis vectors into superpositions: |0) — > 
(|0) + |1))/V2, |1) ^ (|0) - |1))/V2. This gate is of 
use to prepare a specific initial state: when applied to 
every qubit of the system in the ground state |0. . .0), 
it provides an equally-weighted superposition of all basis 
states: 



H . 



H|0...0) = 



1 



2N/2 



E Ml 

...<i„=0,l 



(B4) 



The terms in the sum can be viewed as binary represen- 
tations of all integers from up to 2^ — 1. Thus the state 
(B4) is a superposition of all these integers. When used 
as input of a quantum algorithm, it represents 2^ clas- 
sical inputs. Due to linearity of quantum time evolution 
they are processed simultaneously, and also the output 
is a superposition of 2^ classical results. This quantum 



^'^VNOT is not uniquely defined: the matrices i ^^"^UxXc 
n-n /2) with n = 1, 3, 5, 7 produce NOT, when squared. 



parallelism is a fundamental property of quantum com- 
putation and is responsible for the exponential speed-up 
of certain quantum algorithms. 

Among the two-qubit gates an important one is the 
exclusive-OR (XOR) or controUed-NOT (CNOT) gate: 



CNOT 



/ 1 \ 

10 

1 

V 1 y 



(B5) 



When applied to classical (basis) states it flips the second 
bit only if the first bit is 1. It was shown by Barenco et al. 
(1995) that the CNOT gate together with single-bit op- 
erations forms a universal set, sufficient for any quantum 
computation. In other words, any unitary transformation 
of a many-qubit system can be decomposed into single- 
bit gates and CNOT-gates. This explains the importance 
of CNOT in the quantum information-theoretical litera- 
ture. However, it should be pointed out that almost any 
two-qubit gate (an exception is the classical SWAP-gate) , 
when combined with single-bit operations, forms a uni- 
versal set. 

Let us also mention another useful two-bit gate, the 
controlled phase shift: 



i?(0) 



/ 1 

10 

10 

\ e*"^ 



(B6) 



It shifts the phase of the state |1) of the second qubit 
when the first qubit is in the state |1). (We use an un- 
conventional symbol for this operation, which stresses its 
symmetry with respect to the transposition of qubits.) 



Quantum Fourier Transformation 

As an example we discuss the quantum algorithm for 
a discrete Fourier transformation. N qubits allow us to 
represent the integers j = 0, . . . , 2^ — 1 as basis states 
|0) , . . . , |2^ — l). Starting from a superposition of these 
states with amplitudes Cj and applying the combination 
of controlled phase shifts and Hadamard gates shown in 
Fig. p| one obtains an output: 



2"-l 2"-l 
J=0 



(B7) 



fc=0 



where the output amplitudes Ck and input amplitudes Cj 
are related by the discrete Fourier transform: 



Ck 



1 ■^-^ / inikj " 
2iv 2^ exp ( — ri— I c. 



2^ 



(B8) 
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While in classical computation the time needed for the 
Fourier transform grows exponentially with the number 
of bits N, the quantum algorithm in Fig. takes oc N"^ 
steps. The quantum Fourier transformation was devel- 
oped by Shor (1994) and later improved by Coppersmith 
(1994) and Deutsch (see Ekert & Josza 1996). Its ex- 
ponential speed-up compared to classical algorithms is 
crucial for the performance of Shor's (1994) algorithm 
for the factorization of large integers. 
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FIG. 21. A realization of the Fourier transformation for 4 
qubits, i.e. 2^ coefficients. The phase shifts are 4'jk = 7r/2-'^'°. 



Quantum computation and optimization 

The unitary transformations needed for quantum com- 
putation or any simple quantum manipulation should be 
realized in a particular physical system. For this pur- 
pose they should be decomposed into elementary uni- 
tary gates. This decomposition is not unique and should 
be optimized with respect to various parameters, for in- 
stance, the time, the number of steps, or the complexity 
of the manipulations involved. In some physical systems 
the manipulations involve additional dissipation (as com- 
pared to the idle state) which should be optimized as 
well. In this subsection we discuss how certain unitary 
logic gates can be realized in a spin system with model 
Hamiltonian (|A^). The results are also useful for many 
other physical realizations of qubits (see Sections |l|, |l|) 
since many of them have similar Hamiltonians with sim- 
ilar control parameters. 

First we consider the Hadamard gate (p^). It can be 
performed, up to an overa ll p hase factor, as a sequence 
of elementary operations ( |A4[) and (A5): H oc Ux{a — 
n/4:)Uz{(3 = TT/4)Ux{a = 7r/4). However, it can also be 
performed faster by simultaneous switching of and 



H oc exp 



'2 V2 



(B9) 



The CNOT gate in the model system (A_2) can be im- 
plemented by a combination of 2 two-qubit gates U2b 
( [A^ ) and several single-qubit gates (see also Imamoglu 
et al. 1999): 

CNOT cx (I) C/f (-|) [7| i-n) U,, (-|) x 



Similarly, the controlled phase shift gate ( |B6| ) is produced 
by the sequence 



RW'xu![^^)u2b {-l)u'x 



(Bll) 



One can see from these examples that it takes quite 
a number of elementary gates to perform the CNOT 
or R{(j)) and further optimization is desired. In many 
realizations two-qubit elementary gates are more costly 
(complicated or longer) than single-qubit gates. With 



this taken into account one can ask whether Eqs. (BIC) 



and (Bll) can be reduced to just one two-bit gate? An 



analysis of this problem (Makhlin 2000) shows that it 
is not possible if the two-bit elementary gate (|A6|), pro- 
duced by the ATy-coupling Hamiltonian, is used. The 
result is the same for the Heisenberg spin coupling, but 
the Ising-type coupling, oc UzCFz, allows to achieve R{(j)) 
[and hence CNOT = \^'^ R{-k)H'^] with only one two-bit 
operation: R{cj)) oc C/j(-0/2)C/2(-(/)/2) exp(i(/)crlcr2/4). 
The optimization of certain quantum logic circuits was 
also discussed in connection with qubits based on spins 
in quantum dots (Burkard et al. 1999). 



APPENDIX C: CHARGING ENERGY OF A QUBIT COUPLED 
TO A SET 

The charging energy of the system shown in Fig. |l^ is 
a quadratic function of the charges n and N: 



Hc(n, N,Vn,VN)^ 4Sc,qbn2 + 
+E;^tN{2n - 1) -f 2enVn + cNVn ■ 



SET 



(CI) 



The form of the mixed term oc N{2n — 1) is chosen 
for later convenience. The charging energy scales i?c.qb, 
£^c,SET and Ei^t are set by the capacitances between all 
the islands. Elementary electrostatics shows that they 
can be written as 

^c,qb = e2(Cf T + c\nt + 2Ct)/2A « eV2Cj , 
EcsET = e2(Cf + Cint + Cj)/2^ « eV(4CT) , 

Sint = e^CintM « e2Cint/(2CjCT) . (C2) 
Here we introduced 

A = (Cf + C7j)(Cf T + a„t + 2Ct) 

+ Cint (Cf ^ + 2Ct) « 2CjCt . (C3) 



For simplicity we assumed that the two tunnel junctions 
of the SET have equal capacitances Ct , and the approx- 
imate results refer to the limit C*, Cint, C^^"^ < Ct < 
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Cj, which we consider useful. The effective voltages Vn 
and Vn depend in general on the gate voltages Vg'^'', 
V^^^""" and the transport voltages applied to the SET's 
electrodes. However, for a symmetric setup (equal junc- 
tion capacitances) and symmetrically distributed trans- 
port bias (as shown in Fig. 14), Vn and Vn are controlled 
only by the two gate voltages: 



Vm^V 



SET 



s A 



Vn=V 



SET 



A e ' 



A 



(C4) 



The total charging energy can thus, up to a nonessen- 
tial constant, be presented as a sum of the contributions 
of the qubit (^ , the SE T (O ), and the interaction 



term EintN{2n-l) (cf. Eq.(|5j)). The effect of the SET 
is to renormalize the parameters of the qubit Hamilto- 
nian ( ^.iD : Ec and rig should be substituted by -Ec,qb 
and ng,qb = -eVn/4Ec,qh- 



APPENDIX D: DERIVATION OF THE MASTER EQUATION 

We briefly review the rules for the evaluation of dia- 
grams; for more details including the discussion of higher 
order diagrams we refer to SchocUer & Schon's (1994) 
paper. Typical diagrams, which are analyzed below, are 
displayed in Figs. 02 and 23, The horizontal lines, dis- 



cussed in detail below, describe the time evolution of the 
system governed by the zeroth order Hamiltonian Tip. 
The directed dashed lines stand for tunneling processes, 
in the example considered the tunneling takes place in 
the left junction. According to the rules, the dashed lines 
contribute the following factor to the self-energy S: 



ah 



TTkBTV exp[±|AiL(t-t')] 



sinh^ {t-t'± iS)] 



(Dl) 



where aL = ^/(47r^e^i?T,L) is the dimensionless tunnel- 
ing conductance, /^l is the electro-chemical potential of 
the left lead, and is the high-frequency cut-off, which 
is at most of order of the Fermi energy. The sign of the in- 
finitesimal term i6 depends on the direction of the dashed 
line in time. It is negative if the direction of the line with 
respect to the Keldysh contour coincides with its direc- 
tion with respect to absolute time (from left to right), and 
positive otherwise. For example, the right part of Fig. ^ 
should carry a minus sign, while the left part carries a 
plus sign. Furthermore, the sign in front of i^hit — t') is 
negative (positive), if the line goes forward (backward) 



with respect to absolute time. Finally, the first order di- 
agrams are multiplied by (—1) if the dashed line connects 
two points on different branches of the Keldysh contour. 



N,m,j 
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FIG. 22. Example of a self-energy diagram for an 'in'-rate. 
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FIG. 23. Example of a self-energy diagram for an 'out'-rate. 

The horizontal lines describe the time evolution of the 
system between tunneling processes. For an isolated cen- 
tral island they turn into exponential factors e*^^^*"* \ 
depending on the charging energy of the system. How- 
ever, in the present case the island is coupled to the qubit, 
and we need to account for the nontrivial time evolution 
of the latter. For instance, the upper line in the left part 
of Fig. ^ corresponds to (iV, j|e~^^"(*~* )|iV,j'), while 
the lower line corresponds to (A^-|-l, z'le^^"*^*"* ^\N+\,i). 

In the present problem we assume that the tunneling 
conductance of the SET is low compared to the quantum 
conductance. In this case lowest order perturbation the- 
ory in the single-electron tunneling, describing sequential 
tunneling processes, is sufficient. The diagrams for E can 
be split into two classes, depending on whether they pro- 
vide expressions for off-diagonal {N' ^ N) or diagonal 
[N' — N) elements of E in N . In analogy to the scat- 
tering integrals in the Boltzmann equation these can be 
labeled 'in' and 'out' terms, in the sense that they de- 
scribe the increase or decrease of a given element p*' 
of the density matrix due to transitions from or to other 
iV-states. Examples for the 'in' and 'out' terms are shown 
in Figs. ^ and ^ respectively. 

We now are ready to evaluate the rates in Figs. ^ and 
For example, the 'in' tunneling process in the left 
jimction is expressed as 



exp 



i ( T^N 



sinh 



{At + iS)] 



(D2) 
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exp 



where E^^^ 



[-Ec,set(-^^i 
[Ec,se.t{N2 - A^s.set)^ - Ml N- 
gain for the tunnehng in the left junction in the absence 
of the qubit, and the operators 



A^g,SET)^ - ML A^l] - 

2 1 is the Coulomb energy 



W^l = Wj(iVi) Ho(N2) 



(D3) 



provide corrections to the energy gain sensitive to the 
qubit's state. Here 'Ho{N) is the A^-th block of the Hamil- 
tonian Tictri + Hint (note that Tictri and Hint are block- 
diagonal with respect to N). The indices and 
relate to the left and right side of the tensor product in 
( p^ ) correspondingly. 

The form of the master equation (5^) suggests the use 
of the Laplace transform, after which the last term in 
(5.8) becomes S](s)/5(s). We Laplace transform (D2) in 



N+l< 



< E^^^, i.e. 



we assume 



the regime hs, \Wl^^^\, \Wj^ 
the density matrix p to change slowly on the time scale 
given by h/E^_^_^. This assumption should be verified 
later for self-consistency. The inequalities also mean that 
we choose the operation regime of the SET far enough 
from the Coulomb threshold. Therefore the tunneling is 
either energetically allowed for both states of the qubit 
or it is blocked for both of them. At low temperatures 
{kBT < E^+i) and for E^^^S < ?i we obtain 



^N,m,i'^N+l,m,i I'^J 



n 7=1 JV 



-aLD{E^^,) 



2s- 



iw, 



N 

N+1 



N 

N+1 



(D4) 



where D{E^_^^^) 
Euler's constant 



i 1 + 7 + ln{E^^j^6/h) and 7 w 0.58 is 
The first term of ( p4D is the standard 
Golden rule tunneling rate corresponding to the so-called 
orthodox theory of single-electron tunneling (Averin & 
Likharev 1991). The rate depends strongly on the charg- 
ing energy difference, E^_^^, before and after the pro- 
cess, which in the present problem is modified accord- 
ing to the quantum state of the qubit (the VF-terms). 
At finite temperatures the step-function is replaced by 
e{E) ^ [1 - cxp {-E/kBT)]-^. We denote the full ma- 
trix of such rates F. As has already been mentioned, we 
concentrate on the regime when the leading tunneling 
process in the SET is the sequential tunneling, involving 
only two adjacent charge states, say, iV = and = 1 
(to avoid confusion with the states of the qubit we keep 
using the notation N and + 1). Let us, for example, 
calculate a submatrix of T which originates from the first 



term on the RHS of (D4) and corresponds to the tun- 
neling process A^^A^-f-1, m—>m in the left junction. 
This submatrix is a super-operator which acts on a 
2 X 2-matrix p as 



KTlP = 2'KaLE^+iP + TTOfL [<57Yint, p] 



+ J 



(D5) 



where <5Hint = UoiN + 1) - n^{N) = Ei^ta,. 

The last, logarithmically diverging term of Eq. (D4) 
produces the commutator term in the RHS of the master 
equation (5.S). These terms turn out to be unimportant 
in the first order of the perturbation theory. Indeed, for 
the left junction we obtain the following contribution to 
the RHS of (PI): 



(D6) 



where Hq 



-^{Hq{N) + Ho{N + 1)) and Z^l is a ma- 



trix in N and m spaces. The eigenvalues of the matrix 
Dh are at most of order D{E^^-j^). Neglecting terms 
of order ai^D{E^_^_-^^)Eint in Eq. ([P^), we can replace 
TYo by Hq. Our analysis shows that these neglected 
"coherent-like" terms do not change the results as long as 
q;l I ln{E^_^-^^6/%)\ ^ 1. A similar analysis can be carried 
out for the right timnel junction of the SET. 

Now we can transfer all the "coherent-like" terms into 
the LHS of the master equation. 



1 - ckl-Dl - aRDn 



(dm 
{ dt 



1 



r/5(0, 

(D7) 



and multiply Eq. (D7) from the left by (1 — ai^Di^ — 
ttR-DR)^^ ~ (1 + a^Di^ + ajiDB.) so that the corrections 
move back to the RHS. Since F is itself linear in and 
ttR, the corrections are of second order in a (more accu- 
rately, they are small if a\\n{E^^-^^S/fi)\ ^ 1 for both 
junctions). Thus we drop the "coherent" corrections and 
arrive at the final form of the master equation: 



^4[p(t),Ho]4m.) 



(D8) 



We have shown that under the assumption of 
sufficiently slow dynamics of the qubit and SET, 
-Bj,B,(t/g),^int,FL/R « E}^+,,E^+\ the evolution of 
the system reduces to Markovian dynamics as described 
by the master equation (D8). 



APPENDIX E: QUANTUM DOTS AND SPINS 

Ultra-small quantum dots have also been suggested as 
possible realizations of qubits. A quantum dot is, for 
instance, a small area in a two-dimensional electron gas 
(2DEG) where electrons are confined by the surrounding 
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potential walls created by metal structures on top of the 
systems and applied gate voltages. Usually the size of the 
dot is larger than the Fermi wave length of the 2DEG, 
Af « 30 — 80nm. Then, many electrons are confined 
in the dot. To construct qubits one needs smaller dots 
whose size is or order of Ap. Such dots resemble artificial 
atoms since they contain only a small number of electrons 
filling discrete energy levels. The simplest qubit would be 
formed by an artificial 'hydrogen molecule' with one elec- 
tron confined in a double dot. Of course, a system with 
many electron states is sufficient as well, provided the 
level spacing is large enough such that by proper tuning 
only one level in each dot plays a role. In these systems 
the coherent tunneling of single electron charges, con- 
trolled by external gate voltages, has been demonstrated 
(Oosterkamp et al. 1998), and many of the ideas out- 
lined in Sectional for Josephson junction devices apply as 
well. However, there appear several difficulties: one has 
to push the fabrication to the limits (in contrast to the 
superconducting box, where we exploit the macroscopic 
coherence of the superconducting state), and the system 
is likely to be disordered due to variations in the difficult 
fabrication process. These problems may be overcome, 
e.g., by using chemically well-defined cluster molecules or 
other advanced growth techniques, but still at this stage 
there is not much of an activity to be reported exploiting 
the charge degrees of freedom in quantum dots. 

Spin degrees of freedom have also been proposed as 
candidates for qubits in mesoscopic systems (Loss & Di- 
Vincenzo 1998, Kane 1998). The most obvious advan- 
tage of this approach is the fact that the Hilbert^space 
of a spin-1/2 particle is really two-dimensional E3. In 
contrast, in most other proposals the two-dimensional 
Hilbert space of a qubit is a subspace of a much larger 
physical Hilbert space, to which the system can leak out. 
The proposal of Loss & DiVincenzo (1998) deals with 
spin states of electrons in quantum dots. The dot has to 
be small enough and the temperature low enough that 
the orbital degrees of freedom of the electron are fixed, 
and the only relevant degree of freedom is the electron's 
spin. The spin degree of freedom is quite decoupled from 
the voltage fluctuations of the gates. Therefore the de- 
phasing time of the qubit may be quite long (of order of 
microseconds). 

To manipulate individual qubits one has to apply indi- 
vidual (local) magnetic fields to each spin. An alternative 
method is to apply a homogeneous magnetic field and 
shift electrostatically the position of the spin in the dot 
to areas with different g-factors (DiVincenzo et al. 2000). 
In this way a resonance condition can be controlled for 
individual spins. The two-bit coupling is achieved by 



controlling the tunneling between neighboring quantum 
dots. This may be done by changing the gate voltage con- 
figuration, which changes the potential barrier between 
the dots. In this way the overlap of the wave functions 
of the electrons in the dots is controlled and so is the ex- 
change splitting between the singlet and the triplet states 
of the two spins produced by the Coulomb interaction. 
This may be represented as an interaction Hamiltonian of 
the Heisenberg type, Hint = J{t)SiS2- Since the overlap 
integral depends exponentially on the barriers and can 
be made very small, also the interaction strength J(t) 
can be tuned by gate voltages to exponentially low val- 
ues, completely decoupling the qubits. If turned on, the 
Heisenberg interaction allows one t o p erform the 'square 
root of swap' operation (see Eq. A6), which, together 
with the one-bit gates, constitute a universal set. 

Another setup which exploits the spin states was pro- 
posed by Kane (1998). in this case the qubit is the nu- 
clear spin 1/2 of the phosphorus atom placed as a donor 
in the silicon host. The host nuclei have zero spin. The 
nuclear spins are addressed by gating the shallow bound 
states of electrons around the positively charged donor 
ions. If a high enough magnetic field is applied, the elec- 
tronic spins are completely polarized. The hyperfine in- 
teraction makes the nuclear spins feel an additional ef- 
fective magnetic field proportional to the density of the 
electronic wave function at the nucleus. The wave func- 
tions of the shallow states extend thousands of A away 
from the donors' nuclei. Therefore one can manipulate 
these states by gate voltages, changing their wave func- 
tions' shape and, effectively, the local magnetic field for 
individual qubits. If one allows the shallow states from 
two donor atoms to overlap, an effective exchange in- 
teraction between the nuclear spins emerges. The wave 
function overlap may be again controlled by gates. Thus 
a controllable two-bit coupling is achieved. 

A direct measurement of the quantum state of an in- 
dividual spin is a difficult task. Instead, in both propos- 
als of Loss & DiVincenzo (1998) and Kane (1998), it is 
suggested to translate the spin state into a charge state 
which could be measured easier. Both proposals are in- 
teresting from a physics point of view and have many 
advantages, allowing in principle for a large number of 
phase coherent manipulations. On the other hand, both 
proposals operate at the limit of present day technology 
and even the simplest quantum manipulations still have 
to be demonstrated. 
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